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Abstract 



We present a consistent low energy effective field theory framework for parameterizing 
the effects of novel short distance physics in inflation, and their possible observational signa- 
tures in the Cosmic Microwave Background. We consider the class of general homogeneous, 
isotropic initial states for quantum scalar fields in Robertson- Walker (RW) spacetimes, sub- 
ject to the requirement that their ultraviolet behavior be consistent with renormalizabil- 
ity of the covariantly conserved stress tensor which couples to gravity. In the functional 
Schrodinger picture such states are coherent, squeezed, mixed states characterized by a 
Gaussian density matrix. This Gaussian has parameters which approach those of the adi- 
abatic vacuum at large wave number, and evolve in time according to an effective classical 
Hamiltonian. The one complex parameter family of a squeezed states in de Sitter space- 
time does not fall into this UV allowed class, except for the special value of the parameter 
corresponding to the Bunch-Davies state. We determine the finite contributions to the infla- 
tionary power spectrum and stress tensor expectation value of general UV allowed adiabatic 
states, and obtain quantitative limits on the observability and backreaction effects of some 
recently proposed models of short distance modifications of the initial state of inflation. 
For all UV allowed states, the second order adiabatic basis provides a good description of 
particles created in the expanding RW universe. Due to the absence of particle creation for 
the massless, minimally coupled scalar field in de Sitter space, there is no phase decoher- 
ence in the simplest free field inflationary models. We apply adiabatic regularization to the 
renormalization of the decoherence functional in cosmology to corroborate this result. 



I. INTRODUCTION AND OVERVIEW 



Inflationary models were introduced principally to account for the observed large scale 
homogeneity, isotropy and flatness of the universe in a causal way, independently of detailed 
initial conditions Because of the exponential expansion of an initially small causal 
patch, the inflationary de Sitter epoch dominated by the vacuum equation of state p = 
—e suppresses any classical inhomogeneities in the initial conditions by many orders of 
magnitude, and leads to a primordial power spectrum that is both scale invariant and 
featureless. Most inflationary models assume that the quantum fluctuations which lead to 
this scale invariant spectrum originate in the maximally 0(4, 1) symmetric Bunch-Davies 
state of scalar fields in de Sitter spacetime although the possibility that other states 



may play a role was considered by some authors 
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In the last few years there has 



been a renewed interest in the possible effects of different initial states in inflation 
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fueled largely by the speculation that more 
precise observations of the Cosmic Microwave Background (CMB) might make such effects 
observable, thus opening up the possibility of using CMB observations to probe novel short 
distance physics in the very early universe. 

The primary purpose of this paper is to present a consistent low energy effective field 
theory (EFT) framework for parameterizing such short distance and initial state effects in 
cosmological spacetimes. Although the elements of quantum field theory in curved space 
upon which this EFT framework relies have been known for some time, a comprehensive 
treatment of general homogeneous, isotropic initial states in Robertson- Walker spacetimes 
has not been given previously to our knowledge. Such a treatment of general initial states 
requires both covariant and canonical methods, as well as a dictionary to translate between 
them. Establishing the EFT framework, the relationship between the covariant and canon- 
ical approaches, and the form of the state-dependent terms in the covariant stress tensor 
occupies Sees. ITTllI VI of the paper. The paper is designed so that after becoming acquainted 
with the definitions and conventions in Sec. |TIJ the reader may skip the detailed development 
of Sees. Illim VI if desired, and go directly to the applications in later sections, referring back 
to the previous sections for the derivation of the formulae as necessary. Readers interested 
only in particular modifications of the initial state of inflation and their effects on the CMB 
may wish to skip directly to Sec. El 
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The secondary purpose of this paper is to apply the EFT methods developed in Sees. HTl 
II VI to the processes of particle creation and decoherence in semi-classical cosmology. Certain 
problems with the definitions of particle number and the decoherence functional are resolved 
by the same adiabatic methods used to define the general class of UV allowed initial states 
in the EFT approach. Readers interested primarily in these applications may wish to go 
directly to Sees. IVII or I VI II respectively, likewise referring to the earlier sections as needed. 
The remainder of this Introduction contains a general overview of the issues addressed in 
the paper and our approach to them. As a further guide to the content of the paper, the 
concluding Sec. IVIIII contains a point-by-point summary of our main results. 

As we consider initial state modifications of inflation it is perhaps worth emphasizing 
from the outset that any sensitivity of present day observations to initial state effects runs 
counter to some of the original motivations for and attractiveness of inflation. A scale 
invariant spectrum is one of inflation's most generic predictions, precisely because of the 
presumed late time msensitivity to perturbations of the initial state. If there are features in 
the power spectrum of the CMB today which are not erased by the exponential redshift of 
the inflationary epoch and which bear the imprint of new physics at short distance scales, 
then one might ask what prevents short distance physics from affecting other large scale 
properties of the universe, such as its homogeneity, isotropy or flatness. Since it is not clear 
which inflationary model (if any) is correct, fine tuning a specific model to make particular 
modifications observable in the CMB power spectrum results in a diminishing of the overall 
predictive power of inflation. As long as it is possible to accommodate any observable 
features in the power spectrum by appropriately fine tuning the inflationary model, the 
physical origin of these features as true signatures of new high energy physics must remain 
uncertain |2(|. Finally, the remarkable detection of a non-zero cosmological dark energy in 
the universe today [2j|, at a level very different from estimates based on considerations of 
"naturalness" from short distance physics, should caution that present cosmological models 
are as yet far from complete, and the connection between microphysics and macroscopic 
structure in the universe is still to be elucidated. 

Despite these fine tuning and naturalness problems, it is nevertheless true that in almost 
any given inflationary model there can be surviving initial state effects in the primordial 
power spectrum at some level, and the advent of more precise CMB data makes quantifying 
the sensitivity of inflationary models to such initial state effects potentially worthwhile. 



2 



A quantitative treatment of short distance and initial state effects in inflation requires a 
predictive, low energy framework in which such effects can be parameterized and studied with 
a minimum of assumptions about the unknown physics at ultrashort distances. Effective 
field theory (EFT) provides exactly that framework in other contexts, and we assume in 
this paper that EFT methods may be applied in gravity and cosmology as well. The EFT 
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approach to perturbative gravitational scattering amplitudes was discussed in Ref. 
we extend EFT methods to the non-perturbative regime of semi-classical cosmology. 

Since a fully predictive quantum theory of gravity is still lacking, we are virtually com- 
pelled to adopt an EFT approach to cosmology. Since all scales are presumed to be redshifted 
to lower energy scales where an EFT description eventually becomes applicable, the EFT 
appropriate for cosmology which respects general coordinate invariance and the Equivalence 
Principle is the Einstein theory together with its quantum corrections. In the EFT frame- 
work of semi-classical gravity we can study general (i.e., scale non-invariant) initial states 
in Robert son- Walker (RW) spacetimes in a well-defined way, without detailed knowledge of 
the short distance physics which may have generated them. Although more general initial 
states and more general matter EFT interactions could be considered, we focus in this pa- 
per on the specific gravitational effects of quantum matter fields, and restrict ourselves for 
simplicity to free scalar fields in spatially homogeneous and isotropic initial states consistent 
with the symmetry of the RW geometry. 

The red-shifting of short distance scales to larger ones as the universe expands distin- 
guishes semi-classical gravity from other effective field theories, which possess a fixed physical 
cutoff. In the case of a fixed cutoff, the shorter distance modes of the effective theory can be 
excluded from consideration, and their effects subsumed into a finite number of parameters 
of the low energy description. In practice, absorption of the cutoff dependence of quantum 
corrections into a finite set of parameters of the effective action is no different in an EFT 
from that in a renormalizable theory, except for the allowance of higher dimensional inter- 
actions and the corresponding parameters which are suppressed by the cutoff scale in the 
EFT description. 

Implicit in the EFT framework is the assumption that the effects of short distance de- 
grees of freedom decouple from the long distance ones. However, in an expanding universe 
decoupling is a delicate matter. New short distance modes are continually coming within 
the purview of the low energy description, and some additional information is required to 
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handle these short distance degrees of freedom as they newly appear. If these ultraviolet 
(UV) modes carry energy-momentum, as they should when their wavelength becomes larger 
than the short distance cutoff and EFT methods apply, then simply excluding their con- 
tribution to the energy-momentum tensor at earlier times will lead to violations of energy 
conservation, a well-known point that has been emphasized anew in Ref. ^j]. Energy non- 
conservation occurs with a physical momentum cutoff because energy is being supplied (by 
an unspecified external mechanism) to the new degrees of freedom as they redshift below the 
cutoff, although they carried no energy-momentum formerly. This is an essential point: arbi- 
trary short distance modifications that violate energy conservation are unacceptable in a low 
energy EFT respecting general coordinate invariance, since the resulting energy-momentum 
tensor (T^) cannot be a consistent source for the semi-classical Einstein equations at large 
distances. Upsetting the macroscopic energy conservation law by a short distance physical 
cutoff affects the cosmological evolution at all scales, hence violating decoupling as well. 

These considerations show that the necessary existence of a conserved source for the 
semi-classical Einstein equations provides an important constraint on the class of possible 
short distance modifications of the initial state of inflation, quite independent of the matter 
field content and its EFT. Let us emphasize that there is no problem modifying the initial 
state of a quantum field at a fixed time to for momenta below some physical scale M at 
that time. However the quantum state is not completely specified and a conserved energy- 
momentum tensor cannot be computed unambiguously, until information is given also for 
physical momenta initially much greater than M, which will redshift below M at later times 
t > to. Without some prescription consistent with general coordinate invariance for dealing 
with these arbitrarily high energy "trans-Planckian" modes, which will become physical 
low energy modes eventually, the low energy effective theory of semi-classical gravity is not 
complete or predictive. 

General covariance of the low energy effective theory of gravity coupled to quantum 
matter is the key technical assumption which we make in this paper. General coordinate 
invariance determines the form of the effective action, and therefore the counterterms which 
are available to absorb the ultraviolet divergences of the energy-momentum tensor of the 
quantum matter fluctuations. The renormalization of (T &) with the standard local covariant 
counterterms up to dimension four is possible if and only if the short distance properties of 
the vacuum fluctuations are the standard ones, as expressed for example in the Hadamard 
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conditions on the two-point function as x — > x' |23l.l24j. These UV conditions on 

the structure of the vacuum should be viewed as an extension of the Equivalence Principle 
to semi-classical gravity, since they amount to assuming that the local behavior of quantum 
matter fluctuations are determined in a curved spacetime by those of flat spacetime and 
small (calculable) deviations therefrom. With this physical assumption about the local 
properties of the vacuum, the ultrashort distance modes are necessarily adiabatic vacuum 
modes and are fully specified as they redshift below the UV cutoff scale. Then covariant 
energy conservation is ensured, there are no state dependent divergences in (T ab ), and the 
effective field theory of semi-classical gravity applied to cosmology becomes well-defined and 
predictive within its domain of validity. 

Several authors have considered specific short distance modifications, such as modified 



dispersion relations for the modes 



llj or spacetime non-commutativity 



121. We do not 



consider in this paper these or any other possible specific short distance modifications that 
would take us out of the framework of the covariant low energy EFT of gravity, without 
providing a completely consistent quantum alternative. Once the low energy EFT of semi- 
classical gravity applies, any imprint of UV physics can be encoded only in the parameters 
of the initial state up to some large but finite energy scale M. While many papers have 
discussed possible imprints of new short distance physics on the CMB power spectrum, a 
few authors have considered also the constraints that may arise from the energy-momentum 
tensor of the fields in a state other than the BD state, making use of various order of 
magnitude estimates RU U, HHQ. If the initia! state ^cations ate patametet- 
ized by adding irrelevant higher dimensional operators of the scalar EFT at the boundary, 
there is apparent disagreement between several authors about the order of the corrections 



these modifications induce in the energy-momentum tensor |16l Il7|. Quantitative control of 
the finite state dependent terms in the energy-momentum tensor is potentially important 
for determining whether the short distance modifications can be observable in the power 
spectrum without upsetting other features of inflation. If the scalar field ener gy d ensity is 
too large, it could prevent an inflationary phase from occurring at all In this 

paper we present the framework necessary for the unambiguous evaluation of initial state 
effects in both the power spectrum and the renormalized (T a b) for any homogeneous and 
isotropic state. We then use the general framework to find specific constraints on initial 
states, such as those proposed in the boundary action formalism of Refs. j^j]. We find by 
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2J. Although it had 



explicit computation that the terms in the matter EFT boundary action linear in the first 
higher dimensional operators that one can add do enter the energy-momentum tensor and 
may place more restrictive bounds on the parameters than reported in 

A systematic study of more general initial states in RW spacetimes and their stress 
tensor expectation values was initiated in two recent papers 
generally been assumed that initial state inhomogeneities in (Tab), different from the BD 
expectation value, would redshift on an expansion time scale H~ x , the proper treatment of 
the initially ultrahigh frequency (trans-Planckian) modes is critical to proving this result, 
and also for demonstrating how it may break down in certain special cases. The conditions 
on homogeneous and isotropic initial states of a free scalar field in de Sitter space, necessary 
for its two-point function and energy-momentum tensor to be both finite in the IR and 
renormalizable in the UV were defined in Ref. 2?|. When these conditions are satisfied 
it was shown that (T a b) for a free scalar field of mass m and curvature coupling £ does 
approach the BD value with corrections that decay eu for Rev < 3/2, where ads is 

the RW scale factor in de Sitter space, given by Eq. (|5.1jl and v is defined by Eq. (15 .2j) below. 
For sufficiently massive fields v is pure imaginary and these fields have energy-momentum 
tensor expectation values which decay to the BD value as a~^, just as would be expected 
for classical non-relativistic dust with negligible pressure. However, the light or massless 
cases in which v 2 > 0, show a quite different late time behavior. The massless conformally 
coupled scalar field, for which v = ±1/2, has in addition to the expected subdominant 
a~fg behavior of classical massless radiation with p = e/3, a much more slowly falling a~["g 
component, with p = —e/3, arising from quantum squeezed state effects. In the massless, 
minimally coupled case, relevant for the slowly rolling inflaton field, as well as the graviton 
itself, v — 3/2, and there is an additional a^ s constant coherent component in the late time 
behavior, with p = —e, signaling the breakdown of the 0(4, 1) invariant BD state in the IR, 
and the change in the stress tensor from the BD to the Allen- Folacci ( AF) value j3| . These 
features could not be so readily anticipated by purely classical considerations, but are quite 
straightforward to obtain with the general properly renormalized (Tab). 

We present here a comprehensive treatment of general initial states in arbitrary RW 
spacetimes begun in Refs. HQ, 

from both a canonical and covariant viewpoint. We 
emphasize throughout the paper that in a general RW spacetime all such states are on an 
equal footing a priori. There is no need to resolve the vacuum "ambiguity" often said to exist 
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in curved space field theory. Different physical initial data will simply select different physical 
states. Only the local behavior of these states at very short distances needs to be restricted 
by the Equivalence Principle. Because the structure of the vacuum is most explicit in its wave 
functional representation, we review the Schrodinger description of arbitrary RW states first. 
The Schrodinger description in RW spacetimes has previously been investigated for pure 



states in 



29 



30| and for mixed states in 



] . The general state of a free scalar field is that 
of a mixed, squeezed state Gaussian density matrix p, given by Eq. ()3.10j) . evolving according 
to the quantum Liouville equation ()3.29j) in the Schrodinger representation. The parameters 
that specify this general Gaussian density matrix are in one-to-one correspondence with the 
amplitudes that define the two-point Wightman function and power spectrum of the field. 
We show that the quantum Liouville equation or the scalar wave equation in the covariant 
description imply that these Gaussian parameters evolve with time according to an effective 
classical Hamiltonian (|A10|) . in which h appears as a parameter. This demonstrates that the 
evolution of an arbitrary initial state is completely unitary and time reversible in any RW 
spacetime. Any apparent discrepancy between the Hamiltonian and covariant approaches 
is resolved by including the RW scale factor on an equal footing with the matter field(s) 
in the Hamiltonian description. In any case the correct source for Einstein's equations and 
backreaction considerations is not the canonical or effective Hamiltonian of the Schrodinger 
representation but the expectation value of the covariant energy- momentum tensor, (T a b). 

The fourth order adiabaticity condition on the short distance components of the 
wave functional defines the class of UV allowed RW states consistent with the low energy 
effective field theory satisfying general covariance. We exhibit the finite state-dependent 
contributions to (T ab ) for a general homogeneous, isotropic RW state in (j4.16|) . The 0(4, 1) 
invariant Bunch-Davies (BD) state is a UV allowed fourth order adiabatic vacuum state in 
de Sitter space, but the one complex parameter family of squeezed state generalizations of 
the BD state 0, 0] (sometimes called a vacuua) are not UV allowed RW states . All 
such states except the BD state are therefore excluded as possible modified initial states 
in the low energy description, unless they are cut off at some physical momentum scale 
M j^, and are thence no longer de Sitter invariant. Various possible modifications of 
the inflaton initial state up to some physical scale M at the initial time to are considered 
in Sec. |Vj and their power spectrum and backreaction effects are computed in a consistent 
way. We compare our treatment of short distance and initial state effects with previous work 
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involving a vacuua j8j], adiabatic states jjj, and a boundary action approach 
precise connection between the boundary action approach and the initial state specification 
is established. Readers interested in only these initial state effects in inflation may wish to 
skip directly to Section |3 where these short distance effects are considered, and bounds on 
the short distance modifications are obtained from backreaction considerations. 

In addition to the power spectrum and energy-momentum source for the gravitational 
field, the adiabatic method provides a consistent framework to discuss particle creation 
and decoherence in semi-classical cosmology. Although the definition of particle number in 
an expanding universe is inherently non-unique, we show that the adiabatic number basis 
matched to the second adiabatic terms in (T &) is the minimum one that allows for a finite to- 
tal particle number with conserved energy-momentum. Matching the particle basis to fourth 
or higher adiabatic order is possible but of decreasing physical and practical importance. 
The EFT description also suggests that one should limit the particle number definition to 
the minimal one that requires the fewest number of derivatives of the scale factor, i.e., two, 
which are sufficient to eliminate all power law cutoff dependences in the stress tensor. Hence 
the second order adiabatic basis is selected by the short distance covariance properties of 
the vacuum, together with the local derivative expansion characteristic of a low energy EFT 
description of Einstein's equations, which are themselves second order in derivatives of the 
metric. The second order definition of adiabatic particle number, matched to the form of 
the conserved (T a b) is also the minimal one needed to render the total number of created 
particles in Eq. (jfj.lfjj) finite. As is well known, in the case of the massless, conformally 
coupled scalar field in an arbitrary RW spacetime the zeroth order adiabatic vacuum modes 
become exact solutions of the wave equation, and hence no mixing between positive and 
negative frequency modes occurs. We show that no particle creation occurs also for the 
massless, minimally coupled scalar field (which sometimes serves as the inflaton field) in the 
special case of de Sitter spacetime. 

Particle creation may be described as a squeezing of the density matrix parameters and 
corresponds to a basis in which the off-diagonal elements of p are rapidly oscillating in 
phase, and may be replaced by zero with a high degree of accuracy. To the extent that 
this approximation is valid and the information contained in these rapidly oscillating phases 
cannot be recovered, the evolution is effectively dissipative at a macroscopic level, despite 
being microscopically time reversible. The macroscopic irreversibility is measured by the 
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von Neumann entropy f)6.10j) of the phase averaged density matrix in the adiabatic particle 
basis. In the two special massless cases of cosmological interest mentioned above, namely 
the conformally coupled, scalar field in arbitrary RW spacetime and the minimally coupled 
scalar field in de Sitter spacetime, this phase averaging effect is absent, since no particle 
creation occurs. 

The decoherence functional is defined in the Schrodinger representation as the wave func- 
tion overlap between two states with similar initial conditions but different macroscopic RW 
scale factors, i.e., it measures the quantum (de)coherence between different semi-classical 
realizations of the universe Q]. Physical expectations of a very nearly classical universe 
suggest that this quantum overlap between different macroscopic states in cosmology should 
be finite in principle but extremely small. However, a naive computation of the decoherence 
functional in RW cosmology is plagued by divergences, qualitatively similar to those encoun- 
tered in (T a h) ■ Moreover, previous authors have found that the exact form and degree of these 

nnn 

divergences depend upon the parameterization used for the matter field variables 34, 35, 36J. 
These divergences and ambiguities have prevented up until now the straightforward appli- 
cation and physical interpretation of the decoherence functional in semi-classical cosmology. 
By analyzing the general form of the divergences in the decoherence functional and relating 
them to the divergences in the effective closed time path (CTP) action of semi-classical grav- 
ity (3?!], we show that these divergences can be regulated and removed by a slightly modified 
form of the adiabatic subtraction procedure J3| used to define both the renormalized (Tab) 
and the finite particle number basis. This gives an unambiguous definition of a physical, 
UV finite decoherence functional for RW spacetimes which is displayed in Eq. ()7.13j) and 
which is free of field parameterization dependence and other ambiguities previously noted 
in the literature. The finite decoherence functional does fall rapidly to zero with time in the 
general case in which particle creation takes place, in accordance with physical intuition. In 
the special massless cases in which no particle creation occurs, the renormalized decoherence 
functional vanishes, showing that no decoherence of quantum fluctuations between different 
semi-classical RW universes occurs in these cases. This corroborates the close connection 
between the particle creation and decoherence effects which has been found in other con- 

n 

texts [39J, and shows that the emergence of a classical universe from initial conditions on a 
massless field must be due to other effects, such as interactions, which are neglected in the 
free field treatment presented in this paper. 
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The outline of the paper is as follows. In the next section we establish notation and define 
the general class of homogeneous, isotropic RW states in a RW spacetime. In section IIII1 
we review the Hamiltonian description of the evolution of these states and give the form 
of the mixed state Gaussian density matrix of the Schrodinger representation, as well as 
the Wigner function and effective classical Hamiltonian which describes the evolution. In 
section HVl we evaluate the expectation value of the energy- momentum tensor, and the low 
energy effective action for gravity of which it is part. We specify the conditions on the short 
distance components of a general RW state in order for (T a f,) to be UV renormalizable with 
geometric counterterms of the same form as the effective action, and obtain an expression 
for the finite contributions of arbitrary UV allowed RW states. In section |V] we consider 
three types of modified initial states in inflation, evaluating the power spectrum and energy- 
momentum tensor for each in turn. In section IVII we define the adiabatic particle number 
basis and show how particle creation leads to an effective dissipation in the density matrix 
description. In section I VI II we define a finite renormalized decoherence functional for semi- 
classical cosmology, and corroborate the non-decoherence of massless inflaton fluctuations 
in de Sitter space. We conclude with a detailed summary and discussion of our results. 
Technical details of the Gaussian parameterization of the density matrix and its properties, 
the evaluation of integrals needed in Sec. IS and the comparison of adiabatic bases used in 
squeezing and decoherence calculations by previous authors are relegated to Appendices A, 
B and C respectively Throughout we set c = 1 and use the metric and curvature conventions 
of MTW lioj. 

II. GENERAL RW INITIAL STATES 

Homogeneous and isotropic RW spacetimes can be described by the line element, 

ds 2 = -dt 2 + a 2 (t) dE 2 = -dt 2 + a 2 (t) 7ii dx i dx j , (2.1) 

with t the comoving (or cosmic) time, and 7 y the metric of the three-dimensional spacelike 
sections E of constant spatial curvature, which may be open, flat, or closed. It is also useful 
to introduce the conformal time coordinate, 

r* dt , , 

'=/ wv (2 ' 2) 
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so that the line element (|2.1|) may be expressed in the alternative form, 

ds 2 = a 2 (7?)(-ay + d£ 2 ), (2.3) 

where a is now viewed as a function of conformal time rj. We take a to have dimensions of 
length with rj dimensionless. The scalar curvature is 



R = Q IH + 2H 2 + A,) , H = -. (2.4) 



The overdot denotes differentiation with respect to t and e = —1, 0, +1 depending on whether 
the spatial sections are open, flat, or closed, respectively. 

A free scalar field with mass m obeys the scalar wave equation, 

(-□ + m 2 + £R) $ = 0, (2.5) 

where □ = g V a V& and £ is the arbitrary dimensionless coupling to the scalar curvature. 
Since the RW three-geometry £ is spatially homogeneous and isotropic, the wave equa- 
tion (|2.5)l may be solved by decomposing x) into Fourier modes in the general form, 

$(i,x) = J [dk] (a k fc (t)F k (x) +aj c 0*(t)F k *(x)) . (2.6) 

The Yk are the eigenfunctions of the three-dimensional Laplace-Beltrami operator A 3 on E, 
satisfying 

-A 3 r k (x) ^ V7^ k (x) = (& 2 - e)^k(x) , (2.7) 

and the 4>k{t) are functions only of time and the magnitude of the wave vector, k = |k|. 

For flat spatial sections, e = 0, 7^ = 5^, 7 = det 7^ = 1, and the lk(x) are simply plane 
waves e lkx . The integration measure in Eq. (j2.6J) for this case is J[dk] = J d 3 k/(27r) 3 . 

In the case of compact spatial sections, e = +1, the wave number k takes on discrete 
values which we label by the positive integers k > 1, and the harmonic functions in Eq. (|2.7|) 
are the spherical harmonics of the sphere S 3 . These S 3 harmonics denoted by Yum depend 
on three integers k ^ (k,£,m), the first of which may be identified with |k|, while (£,m) 
refer to the usual spherical harmonics on S 2 with £ < k — 1. Since X^=o (2^+ 1) = k 2 , a given 
eigenvalue of the Laplacian (j2.7J) labelled by k is /c 2 -fold degenerate. The scalar spherical 
harmonics Ym™. may be chosen to satisfy ^(x) = F_ k (x) = F/^_ m (x), and normalized on 



S" 3 so that 



4l| 



/ d 3 E F fe %, m /(x) Y kim (x) = 5 kk , = 5 kk ,5 w 5 mm , , (2.8) 

JS3 
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and 



fc-i 



£=0 m= 



2tt 2 ' 



(2.9) 



which is independent of x. 

In the open case e = — 1, the sums over (£, m) remain, but k becomes a continuous variable 
with range [0, oo). After integration over the direction of k in the e = case, one is also left 
with the integration over the magnitude k with the scalar measure J dk k 2 / (2-7T 2 ). Because of 
Eq. ()2.9|) . the compact e = +1 case is simply related to the non-compact cases of e = 0, —1 by 
the replacement of the integral J dk k 2 / (2n 2 ) by the discrete sum, Ylk=i ^ ' I (2 71 " 2 )- Beginning 
the sum from k = 1 (so that the spatially homogeneous mode on S 3 has eigenvalue k = 1 
instead of k = 0) makes this correspondence between the discrete and continuous cases most 
immediate. We define the scalar measure, 



[d*] = 



/ °°dA; if e = 0,-l 
Er=i if ^ = 1 



(2.10) 



in order to cover all three cases with a single notation. 

The time dependent mode functions 4>k{t) satisfy the ordinary differential equation, 



d 2 fc d0 fc (k 2 - e) 2, fm , n 

d^ dt a" 1 

3 



(2.11) 



If one defines fk(t) = a 2 0fc, then this equation is equivalent to 



fk + 



n „ i / .v 



where 



fc 2 

o; 2 (t) = -+m 2 . 



f k = o 



(2.12) 



(2.13) 



Eq. ()2.12j) is the equation for a harmonic oscillator with time dependent frequency. Note that 
the comoving momentum index k of the mode is constant, while the physical momentum 
p = k/a redshifts as the universe expands. 

An analogous time dependent harmonic oscillator equation may be derived also in con- 
formal time under the substitution, Xki.il) = a 4>k, viz., 



r + mV + (6£ - 1) — + e 



Xk = o. 



(2.14) 



where the primes denote differentiation with respect to conformal time rj. 
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In view of the completeness and orthonormality of the spatial harmonic functions lk(x), 
it is easily verified that the equal time commutation relation 

= ^5 E (x,x') = ^ v ^5 3 (x-x'), (2.15) 
a 6 a 6 

is satisfied, provided the creation and annihilation operators obey 

[ Q k> a k'] = <W > ( 2 - 16 ) 
in the discrete notation, and the complex mode functions satisfy the Wronskian condition, 
a\^ k - ^l) = f k ft - fj* = X ' k xl ~ X k X*k = • (2-17) 

From the equation of motion ([2.11)1 . ()2.12|) or ()2.14|) this Wronskian condition is preserved 
under time evolution. Hence any initial condition for the second order equation of motion 
satisfying (J2.17)) is a priori allowed by the commutation relations. Given any two solutions 
of (|2.11|) . we define their Klein-Gordon inner product as 

io? 

(y, fc ,&) = _ (V^-VM, (2-18) 

which is independent of time. 

Let Vkif) be some particular set of time dependent mode functions satisfying Eq. ()2.11|) 
and the Wronskian condition 1)2.17)1 . These can be used to define a vacuum state. Any other 
set of solutions (f) k satisfying the same Wronskian condition can be expressed as a linear 
superposition of v k and its complex conjugate, 

<P k = A k v k + B k v* k , (2.19) 

which is the form of a Bogoliubov transformation. Because of Eq. (j2.17j) the time indepen- 
dent complex Bogoliubov coefficients must satisfy 

\A k \ 2 -\B k \ 2 =l, (2.20) 

for each k. This is one real condition on the two complex numbers A k and B k . Since 
multiplication of both A k and B k by an overall constant phase has no physical consequences, 
there are only two real parameters needed to specify the mode function for each k. 

The inner product (|2.18J) is preserved under the Bogoliubov transformation (|2.19|) . i.e., 

(<Pk,<frk) = (v k ,v k ) = 1, (2.21a) 
(<f>* k ,<f>k) = (vt,v k ) = 0, (2.21b) 
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$(*,x),— (t,x') 



Thus, we can invert ()2.19|) and solve for the Bogoliubov coefficients at an arbitrary initial 
time, t = t or r\ = r] , with the result, 

A k = (v k , (f> k ) = -£(v* k 0, - vt 0,) o , (2.22a) 

B k = (v* k , 4> k ) = — Mv k <l> k - v k <f) k ) . (2.22b) 
n 



Interactions may be incorporated in this treatment as well |42t |44j . within the semi-classical 
large N approximation, but in order to keep the discussion as simple as possible we shall 
not consider scalar self-interactions in this paper. 

We restrict our attention to initial states of a free scalar field, which like the RW ge- 
ometry (j2.1j) itself, are also spatially homogeneous and isotropic, and call such states RW 
states. Spatial homogeneity of the RW states, i.e., invariance under spatial translations in 
E implies that the bilinear expectation values, (a k a k /) and (a k a k /), can be non- vanishing if 
and only if k = k', while the expectation values, (a^a^i) and (a^a k i), can be non-vanishing if 
and only if k = — k'. In addition, isotropy of the RW states under spatial rotations implies 
that the expectation value of the number operator for k = k', 

( a W) =n k = (a k al) - 1 , (2.23) 

can be a function only of the magnitude k. This constant number in each Fourier mode is 
the consequence of the unmeasurable U(l) phase of the mode functions, and becomes the 
third real parameter needed for each k to specify the initial quantum state of the scalar field. 
We show in the next section that if the free field density matrix for a RW state is described 
by the general Gaussian ansatz in the Hamiltonian description, then the state is necessarily 
a mixed state if n k 0. 

Because of the two parameter freedom to redefine 4> k according to the Bogoliubov trans- 
formation, 1)2.19)1 and (J2.20D . it is always possible to fix the parameters so that the remaining 
bilinears are equal to zero [44j . i.e., 

( Gk a_ k ) = (4a t _ k ) = 0. (2.24) 

If $ is expanded in terms of the vacuum modes v k instead of <p k , then the corresponding 
annihilation and creation operators are 

c k = A k a k + B* k al k , (2.25a) 
4 = A* k a{ + B k a_ k . (2.25b) 
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This is the characteristic form of a Bogoliubov transformation to a squeezed state. If the 
arbitrary overall phase is fixed by requiring A k to be real, then the general squeezed state 
parameters r k and 6 k are defined by 

A k = coshr fc , (2.26a) 
B k = e idk sinhr fc . (2.26b) 

The bilinear expectation values, 

(c k c_ k ) = (2n k + 1) A k B* k = a k A k B* k , (2.27a) 

(c k cl k ) = (2n k + 1) A* k B k = a k A* k B k , (2.27b) 

are non-zero, and 

Nk = (4c k ) = \B k \ 2 + (2\B k \ 2 + l)n k = n k + a k \B k \ 2 , (2.28) 

is the average occupation number of the general mixed, squeezed state with respect to the 
vacuum modes v k . We have introduced the shorthand notation a k = 2n k + 1 for the Bose- 
Einstein factor in Eqs. fTI7\ and (J22BJ)- 

The two-point (Wightman) function of the scalar field may be expressed in terms of the 
mode functions <p k and n k in the form, 

($(t,x)$(t',x')> = J [dk] rn fc 0;(O^(O + (n* + l)^(O«(O)>kW>k(>O. (2-29) 

in the case that the expectation value (<&(£, x)) = 0. When x = x', the sums or integrals 
over the angular part of k can be evaluated with the result that 



is independent of x. Here 



($ 2 (t,x)>= / [d*]5M, (2.30) 

P^k;t) = ^a k \Mt)\ 2 (2.31) 

is the power spectrum of scalar fluctuations in the general RW mixed, squeezed initial state. 
It may also be expressed in terms of v k in the form, 

Pt(k;t) = P v {k ] t) + ^[N k \v k {t)\ 2 + a k Re[A k B* k vl(t)^j (2.32a) 

where 

P v (k;t) = ^\v k (t)\ 2 , (2.32b) 
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is the fluctuation power spectrum in the selected vacuum state, and Eqs. (j2.19j) and (J2.28JI 
have been used. 

In inflationary models this scalar power spectrum becomes the source for scalar metric 
fluctuations, and it is the power spectrum of these metric fluctuations that is actually mea- 
sured in the CMB [53]. Due to the linearity of the metric fluctuations and the assumed 
spatial homogeneity of the classical inflaton field, the resulting power spectrum in the CMB 
is the same (up to an overall normalization) as the quantum scalar field spectrum ()2.32|) 
that generates it. For example, in slow roll inflationary models, the relation between the 
linearized curvature perturbation TZ^ and the quantum scalar field perturbation 50k is jl| 

H „ , 



2.33) 



-i t=t» 



where <p(t) is the classical inflaton field (assumed independent of position) and £* is a time 
a few e-folds after the perturbation has exited the horizon. From (J2.33|) we find that the 
power spectrum of Gaussian curvature fluctuations Pn(k;t*) that is actually observed in 
the temperature fluctuations of the CMB is related to the power spectrum of scalar field 
fluctuations P^k; t*) by 

= J_ (JL) 2 P ^ (2 34) 

~ 87T 2 e \M Pl ) P BD (M*) ' K ] 

In the latter relation we have introduced the standard definition of the slow roll parameter 
e Q], not to be confused with the e = 0, ±1 defined in ()2.4|) denoting flat, closed or open 
spatial sections. We have also normalized the spectrum to the scale invariant Bunch-Davies 
vacuum power spectrum P v {k; t*) = P® D (k; £*), given explicitly by Eq. (j5.fi |) of Sec.lVlbelow. 
Since our intention in this paper is to address the short distance and initial state effects of 
the scalar field in a general, model independent way, we focus on the scalar field power 
spectrum P^k; t) of Eqs. ()2.32|) exclusively in the succeeding sections, leaving the model 
dependent connection to the scalar metric power spectrum Pn{k\ t) unspecified. We do not 
discuss tensor perturbations of the metric at all in this paper. 
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III. DENSITY MATRIX AND HAMILTONIAN 



The development of the previous section in terms of expectation values of the Heisen- 
berg field operator x), with initial data specified by the time dependent complex mode 
functions <pk is well suited to a treatment of the covariant energy-momentum tensor and 
its renormalization. We consider this in detail for general RW states in the next section. 
However, the evolution of a quantum system from given initial conditions may be expressed 
just as well in the Schrodinger representation, and it is the latter approach which makes 
most explicit the specification of the initial state and its unitary evolution in configuration 
space. The Schrodinger representation is also the one best suited to discussions of decoher- 
ence of cosmological perturbations and the quantum to classical transition in inflationary 
models, which we discuss in Sec. I VI II We present the Hamiltonian form of the evolution and 
corresponding density matrix for a general RW initial state in this section, demonstrating 
its full equivalence with the covariant formulation. 

For a free field theory (or an interacting one treated in the semi-classical large N or 
Hartree approximations) it is clear that the two-point function (J2.29j) contains all the non- 
trivial information about the dynamical evolution of the general RW initial state. Even if 
non-zero, higher order connected correlators do not evolve in time either in a free field theory 
or in the leading order large N approximation to an interacting theory. Hence at least as far 
as the time evolution is concerned, a time dependent Gaussian ansatz for the Schrodinger 
wave functional (or density matrix) of the scalar field about its mean value is sufficient in 
both cases. A proof of the equivalence between the large N semi-classical equations and 
the evolution of a Gaussian density matrix has been given for flat Minkowski spacetime 
in Refs. QQ- 

The general RW case differs from the flat spacetime case mainly by the 
appearance of the time dependent scale factor, so that with the appropriate modifications a 
mixed state Gaussian density matrix also exists for the scalar field evolution in cosmology. 

In order to obtain this Gaussian density matrix, let us first derive the Hamiltonian form 
of the evolution Eqs. (|2.5jl and (j2.11j) . In its Hamiltonian form, the variational principle 
both for the matter and metric degrees of freedom, should begin with a classical action that 
contains only first time derivatives of both $ and a. However, the standard classical action 
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for the scalar field plus gravity system, viz., 

S*+S E h = ~ f d^vqitV^g^V^l + m^^^ + ^i- J & 4 x^(R-2k) 

(3.1) 

contains the Ricci scalar R in both the matter action Sq> (when £ 7^ 0) and the Einstein- 
Hilbert action Seh- For RW spacetimes R contains second order time derivatives of the scale 
factor in the H term of Eq. (|2.4jl . In order to remove these second order time derivatives from 
the action one should add a surface term to the action functional above, thereby replacing 
S$ + S E h in Eq. (j3.ip by 

dt I I 8ttG n 



S cl [$,$;a,d] = S$ + S EH + 3 I dtd 3 S " 



(3.2) 



which modifies both the matter and gravitational parts of the classical action at the end- 
points of the time integration, but otherwise leaves the Lagrangian evolution equations away 
from the endpoints unchanged. In fact, it is this classical action S c \ modified by the surface 
term in Eq. (|3.2|) . and not S& + Seh whose Euler-Lagrange variation (which by definition 
has vanishing 5$ and 5a at the endpoints) leads to the scalar field equation of motion (12. 5 j) . 
as well as the Friedman equation for the scale factor. The surface term addition to the 
gravitational action for a general spacetime has been given in Ref. 

With this corrected classical action S c \, the conjugate momentum for the scalar field is 

ILj, = ^ = a 3 ($ + 6£iK>). (3.3) 

If we ignore the Friedman equation for the scale factor for the moment, treating a(t) as an 
externally specified function of time, then the classical Hamiltonian density of the scalar 
field alone is 

n 2 

n$ ee $n $ - £ cl = -4 - 3f #(n*$ + $n$) 

2a 6 

3 

+ ^{d^){d^) + a -m 2 & + 3£a 3 [(6£ - l)H 2 + ^] $ 2 , (3.4) 

where we have symmetrized the second term in this expression involving Il$$, in anticipation 
of the replacement of $ and 11$ by non-commuting quantum operators. 

In the Hamiltonian framework the three independent symmetric quadratic variances, 
($ 2 ), ($11$ + 11$$) and (II|) at coincident times determine the Gaussian density matrix 
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p. The one antisymmetric variance, ($11$ — n<j,$) is fixed by the canonical commutation 
relation, 

[$(t,x),n*(t,x')] =^5 s (x,x'), (3.5) 

which using Eq. ()3.3|) and [$, $] = [11$, H$] = is equivalent to Eq. f)2. 15j) . Let us introduce 
the definitions, 

<y k = 2n k + 1 , (3.6a) 

(k{t) = ^fo~k \4>k\ , (3.6b) 

n k {t) = a 3 (C fc + 6^Cfc) , (3.6c) 

for the time independent Bose-Einstein factor, a k , and the two real functions of time, ( k (t) 
and TT k (t). We show in Appendix A that these definitions allow us to express the three 
bilinear Fourier field mode amplitudes in the form, 

cr k \(pk\ 2 = Cfc , (3.7a) 
a k Re (0fc0 fc ) = ( k Ck , (3.7b) 
-M 2 = C 2 k + ^, (3.7c) 

and this allows in turn for the three independent Gaussian variances at coincident spacetime 
points to be written as 

($ 2 ) = ^| [dk]k\l (3.8a) 

($n$ + n$$> = -L J \Ak] k 2 , (3.8b) 
(nl) = / ^ fc2 (*2 + |f ) • ( 3 - 8c ) 

Thus, the three independent bilinears depend on a set of three real functions of k, (( k , n k ; a k ), 
as expected from our discussion of the initial data in the Heisenberg representation of the 
previous section. The usefulness of this particular set is that (( k , ir k ) will turn out to be 
canonically conjugate variables of the effective Hamiltonian that describes the semi-classical 
time evolution of the general Gaussian density matrix, while a k is strictly a constant of the 
motion. Notice also that the power spectrum defined in (j2.31|) can be written in terms of 
the Gaussian width parameter ( k (t) directly as 

PAk;t) = ^( 2 k (t), (3.9) 
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for both the pure and more general mixed state cases. It is independent of x and the 
direction of k by the spatial homogeneity and isotropy of the RW state. 

The Hamiltonian and corresponding pure state Schrodinger wave functional for scalar 
field evolution in cosmology has been previously given in Ref. Q, |30j]. However, a pure 
state Gaussian ansatz for the wave functional imp oses a constraint on the three parameters 
(Cfc; ^k', Cfc), in fact implying <Jk = 1 for all k ^j. . To remove this restriction one must 
allow for the Gaussian ansatz also to contain mixed terms, so that the density matrix 
p 7^ in general. By simply keeping track of the powers of a(t) and its derivatives, 

it is straightforward to generalize the Minkowski spacetime density matrix to the RW case 
with the result, 



(Q\p(t)W) = (%|/3o(0,p;Co,7r o ;cTo)|go) x n^l^' 7 ^^)'^) 

= po n( 2 <) 4 e -p {-4^ (m 2 + k\ 2 ) 

where the {q^} are the set of complex valued Fourier coordinates of the scalar field which 
are time independent in the Schrodinger representation, i.e., the matrix elements of the 
Heisenberg field operator $(to,x) at an arbitrary initial time to are defined by 

(g|$(t ,x)|g'> = (| [dk]F k (x)g k ) (q\q') . (3.11) 

The latter matrix element is non-vanishing only for g k = g k and is defined precisely by 
Eq. ()3.19|) below. In this Schrodinger coordinate representation the action of the conjugate 
momentum operator 11$ is given by 

(g|n $ (t ,x)|g') = -ih (| [dk]y k (x) A) {q \ q >) . (3.12) 

Since $ is a real field, the complex coordinates are related by g k = g_ k and occur in conjugate 
pairs. Hence we have the rule, 

jrr = ^ (3-13) 

and the ±k terms in the density matrix (|3.1U|) are identical, and may be combined. The 
Qk=o field coordinate is real, and we have separated off the k = component of the density 
matrix p in Eq. (|3.10j) . denoting it by p . In this spatially homogeneous and isotropic Fourier 
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component we may also allow for the possibility of a non-vanishing real mean value of the 
scalar field, 

/oo 
dqoqo (q \p {4>,P; Co, vr ; a )\q ) , (3.14) 
-oo 

which because of the RW symmetry is a function of time only. This spatially homogeneous 
expectation value is the classical inflaton field in inflationary models. Separating this mode 
explicitly from the rest is possible strictly only in a discrete basis, such as that corresponding 
to closed spatial sections, e = +1. The density matrix in the spatially homogeneous sector 
is 

Po(t) = (go|po(0,P;Co,7r o ;o- o )|go) 

= (2vrC 2 H exp{* | (q - q' ) - [(q - 0) 2 + (q' - 0) 2 ] 

+^ i(9b - w - Wo - m + 4?^ - Wo - *)) > ( 3 - 15 ) 

where 

p{t) = a 3 (}(t) + 6fi70(t)) (3.16) 

is the momentum conjugate to the spatially homogeneous mean field 4>(t). 
Real field coordinates (q^, gjj for the k ^ modes may be introduced by 

<fc = ^(<£ - «£) , MO (3.17) 

and the functional integration measure over the field coordinate space defined by 

[-Dq}^dq l[dq^dq{. (3.18) 
loo 

The inner product appearing in (|3.11|) is defined by 

(q\q'} = %o ~ <Q \{ 6(<£ - <£)6{<& - <£) (3.19) 

k>0 

and the general mixed state Gaussian density matrix (|3.10|) is properly normalized, 

Trp= f[Vq)(q\p\q) = l, (3.20) 

with respect to this measure. 
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It is clear from f)3. lUj) that if cr*. = 1 for all k, then the mixed terms vanish and the density 
matrix reduces to the product, 



k 

characteristic of a pure state, with 



. 2 \-I J l<?k| 2 . . |<?k' 2 



(<&) = *(?-k) = (2< fc T* exp <{ -J^L. + iwk ^_ [> , k ^ , (3.22a) 

.2x-i J ■ p / T\ (go - 4>) 2 , . Wo 



*(g ) = (20-4 exp ji | (g - 0) - + ttr j , k = , (3.22b) 

which is the Gaussian pure state Schrodinger wave functional in the Fourier representation 
of the complex field coordinate basis (j3.11|) . The pure state case corresponds to = 0, 
o"fc = 1, and requires only the two real functions of k and t, (£&, n^) for its full specification. 

The Wigner function(al) representation of the Gaussian density matrix is obtained by 
shifting gk — ► <?k + ^k/2 and g k — > q^ — in (j3.10|) . and Fourier transforming p with 

respect to the difference variables Xk, viz., 

F w [q,p] = J [T>x] Yl (2vrn)~ 1 exp ^-^p k x k J («k + y P 4k - y) 

= F (g ,Po) JJ^ftofc^exp (— ^j- - |CfcPk - vr fe g k | 2 ) , (3.23) 
k^o ^ ^* n a k ) 

where p k = p_ k , and 

^o(go,Po) = (7r^(7 )~ 1 exp |- ^° 2 ^ ~ [Co(Po-p) -7ro(go -0)] 2 | (3.24) 

is the Wigner function in the spatially homogeneous k = sector. Note that the normaliza- 
tion of the Gaussian Wigner functional is constant in time, as required for a Hamiltonian 
evolution in phase space. For a given q^ this Gaussian function is peaked on the phase space 
trajectory, 

Pk « ? g k , k ^ , (3.25a) 

Po « P + ^(go-0), k = 0, (3.25b) 
Co 

becoming very sharply peaked on this trajectory in the formal classical limit ft — ► 0, er^ 
fixed, although the width of the peak becomes larger for mixed states with larger Uk (with 
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h fixed) 46] . The Wigner functional (|3.23|) is positive definite for Gaussian states and may 
be interpreted as a normalized probability distribution for any ha k j^- 
The functional integration measure (J3.18|) implies an inner product, 



tfaltfi) = / [Vq\ (* 2 |g>(?|*i> =exp(ir 12 ) 



(3.26) 



between pure states, and a coherence probability functional, 



Tr(pip 2 ) 



W / W] {q\PiW) Wlhlq) = exp (-2T 



12 J 



(3.27) 



for general mixed states in the Schrodinger picture. In the case of pure states the real 
functional T 12 becomes Im T 12 of (J3.26|) . In the case fi\ = p 2 , performing the Gaussian 
integrations in the coordinate representation gives 



Tr(p 2 



[Vq\ / [Vq'\ (q\p\q') (q'\p\q) 



n 



<7k 



exp (~^J [dk]k 2 In a h ^j < 1 



(3.28) 



The inequality is saturated if and only if o k = 1 for all k, in which case p = and 
the equality is simply a consequence of the normalization condition, = 1 on the pure 

state wave functional. If for any k, a k > 1, Tr p 2 < 1, which is characteristic of a mixed 
state density matrix. 

For either a pure or mixed state the Gaussian density matrix satisfies the quantum 
Liouville equation, 

ih^=[H*,p], (3.29) 

provided the time dependent parameters (<j>,p; Ck, K k ) appearing in p satisfy the first order 
equations, 



V 



C 



p = 6^Hp — a 



m 2 + 6£— + 6£(6£ - 1)H 2 



Ti k = 6£Hn k - a 3 



k 2 



+ m 2 + (6£-l) (A, + Q^H : 



K 2 a 2 
4^1 



(3.30a) 
(3.30b) 
(3.30c) 

(3.30d) 
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The first two of these equations are equivalent to the second order equation for the spatially 
homogeneous mean field, 



while the latter two may be combined to yield the second order equation for the Gaussian 
width parameter, 



The last equation is derived in Appendix A from Eq. (j2.11j) and the defining relation (j3.6b|) . 

This establishes the equivalence between the Hamiltonian evolution of the Gaussian den- 
sity matrix ()3.10|) . according to ()3.29|) . and the Lagrangian evolution of general RW initial 
states described in the previous section. Since the Hamiltonian ()3.4j) is Hermitian with 
respect to the field coordinate measure ()3.18|) . the time evolution of p is unitary and the 
normalization ()3.20|) is preserved. Hence there is no dissipation in the system and the evolu- 
tion is fully time reversible in principle. The time evolution of the density matrix parameters 
(4>,p; Cki^k) m ay also be derived from an effective classical Hamiltonian, H e ^ = Tr(7i$p), 
given explicitly by Eq. (JAlOj) of Appendix A. This effective Hamiltonian is just the expecta- 
tion value of the quantum Hamiltonian Tiq, in the general Gaussian state, in which h appears 
as a parameter. Notice that the role of the h term in Eq. ()3.32j) is to act as a "centrifugal 
barrier" for the coordinate (k, preventing the Gaussian width parameter from ever shrinking 
to zero. This width depends on the product ha^, so that the classical high temperature 
limit h<Jk — > ksT/uk is treated on the same footing as the quantum zero temperature limit 
ft/Jk — > h. Both classical thermal and quantum uncertainty principle effects contribute to 
the width of the Gaussian in general. 

The Hamiltonian evolution and the density matrix description of RW states is not mani- 
festly covariant under general coordinate transformations, depending as it does on a partic- 
ular slicing of the four dimensional geometry into three dimensional surfaces S. Since initial 
data must be specified on such a spacelike Cauchy surface, this is the natural 3 + 1 splitting 
for initial value problems in RW cosmology. The equation of motion (|2.5|) is certainly in- 
variant under general coordinate transformations, whereas the initial data must be specified 
on some three surface £ for any particular physical initial state. 

We note that the canonical effective Hamiltonian generating the time evolution of the 
wave functional in the Schrodinger representation is not simply related to the expectation 



+ 3#0 + (m 2 + iR)ct) = , 



(3.31) 




(3.32) 
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value (T u ) of the covariant energy-momentum tensor. The key point in reconciling the 
canonical and covariant energy densities is that the full system of matter plus metric fields 
must be taken into account. After the addition of the surface terms to the standard classical 
action in Eq. (J3.1j) to remove the second derivatives of the metric, the momentum conjugate 
to the scale factor a is 

r q q 

n a = — - = — da + 6£aa$ 2 + 6£a 2 $$ . (3.33) 

oa AttGn 

Then the total Hamiltonian density, constructed in the canonical prescription is 

Wtot = + U a a - £ cl 

' a 3 G tt + a 3 T tt , (3.34) 



8nG N 

where 



G n i{H 2 + -^) . (3.30) 



is the tt component of the Einstein tensor, G a b = Rab — g a bR/2 and 

1 • • 1 TfP 1 / f \ 

T tt = -$ 2 + 6£#$ $ + —YHd&Xdj*) + — $ 2 + 3^ [H 2 + -J $ 2 , (3.36) 
is the component of the covariant energy-momentum tensor 

T ab = -^^^ = (V $)(V 6 $)-^(V c $)(V d $) 

9 

-2^V a ($ V,$) + 2^ cd V d ($V c <l>) + £G a 6<l> 2 - — ^$ 2 . (3.37) 

Invariance under transformations of the time coordinate leads to the classical equation of 
constraint, 

H tot = . (3.38) 

Because of Eq. (13.34)1 this coincides with the tt component of the classical Einstein equations, 
viz., the Friedman equation for RW cosmologies. This constraint equation is equivalent to the 
requirement that the classical theory be invariant under arbitrary time reparameterizations, 
t — > t'(t), a condition which Hamiltonian evolution in a fixed background a(t) does not 
require. The three momentum constraints of spatial coordinate transformations, equivalent 
to the ti components of Einstein's equations are automatically satisfied in any homogeneous, 
isotropic RW state. 
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Hence the Hamiltonian and covariant approaches agree, only for the full system of grav- 
ity plus matter, i.e., provided that the RW scale factor is treated as a dynamical degree 
of freedom, on the same footing as the matter field $. In contrast, the Hamiltonian of 
Eq. ()3.4|) generates the correct evolution of the $ field in a fixed RW background, whether 
or not the scale factor a(t) satisfies Einstein's equations, and with no requirement of invari- 
ance under reparameterizations of time. It is the covariant energy-momentum tensor T a b 
that is conserved and should be used whenever the full cosmological theory of matter and 
gravitational degrees of freedom are under consideration. With this important proviso the 
canonical and covariant formulations of the initial value problem are completely equivalent. 



IV. ENERGY-MOMENTUM TENSOR OF UV ALLOWED RW STATES 

In the previous sections we have defined and described general homogeneous and isotropic 
RW initial states, with no restriction on the set of three density matrix parameters ((k, ^k] &k) 
which describes the state and its evolution. However, because of the Wronskian condi- 
tion ([2. 17)) . that enforces the canonical commutation relations of the quantum field, the 
state parameters do not approach zero rapidly enough at large k for the integrals in ()3.8|) or 
the expectation value of the covariant energy- momentum tensor (Tab), to converge. Hence 
these expressions are purely formal, and a definite renormalization prescription is necessary 
to extract the finite state dependent terms. This is a necessary prerequisite for any dis- 
cussion of short distance, initial state, or backreaction effects in inflation, at least within a 
conventional effective field theory framework. 

Because the energy-momentum tensor is an operator of mass dimension four, it contains 
divergences up to fourth order in the comoving momentum cutoff ku- Requiring that the 
forms of the integrands at large k match those expected for the vacuum up to fourth order 
in derivatives of the metric, allows for all the divergences in (T a b) to be absorbed into coun- 
terterms of the relevant and marginally irrelevant terms of the local gravitational effective 
action (3^. This adiabatic order four condition on the initial state imposes restrictions on 
the set of parameters n k ; cr fc ) at large k, and guarantees that the renormalized expecta- 
tion value (T a b) R will remain finite and well-defined at all subsequent times Conversely, 
failure to impose these short distance restrictions on the initial state leads to cutoff depen- 
dence which cannot be identified with covariant local counterterms up to dimension four in 
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the gravitational action, and which violate the assumptions of a low energy EFT for gravity 
consistent with the symmetries of general covariance implied by the Equivalence Principle. 

The available counterterms up to dimension four in the coordinate invariant effective 
action for gravity are the four local geometric terms A, R, R 2 and C abcd C abcd (the square of 
the Weyl conformal tensor), which can be added to the one-loop action of the scalar field, 
[g] . Hence the low energy gravitational effective action is formally 

S cS [g] = S« [g] + ^J— J d*x ^(R-2A)- 1 -j d 4 x ^ {a C abcd C abcd + (3 R 2 ) , 

(4.1) 



where 



SW [ g ] = y Tr ln(- □ + £R + m 2 ) , (4.2) 



and A, Gn, a, and (3 are bare parameters which are chosen to cancel the corresponding 
divergences in if>W[gr]. A fully covariant renormalization procedure is one that removes all 
divergences in S^fg] by adjustment of the scalar parameters A, Gn, a, and (3 of (|4.1|) . 
and only those parameters, in such a way that the total effective action S^fflj?] and the 
renormalized energy- momentum tensor derived from it, 



is finite (i.e., independent of the cutoff k^) and covariantly conserved: 

^(7^ = 0. (4.4) 

Thus the renormalized expectation value (T a b) R is strictly well defined only by reference to 
the full low energy effective action SW^g] and the equations of motion of the gravitational 
field following from it, 

1 



8nG N 

of which it is a part. 

The local conserved tensors, 



(G ab + Ag ab ) = (T ab ) R + a R ^H ab + (3 R «# a6 , (4.5) 



1 S f 

^H ab J d*x^R 2 = 2g ab nR-2V a V b R + 2RR ab -^R 2 , (4.6a) 

^H ab = j_ / / d 4 ^ C abcd C abcd = AV c V d C acbd + 2R cd C acbd , (4.6b) 

V ~9 °g J 
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derived from the fourth order terms in the effective action are similar to those which appear in 
any EFT, whose equations of motion involve a local expansion in the number of derivatives. 
Provided that we restrict our attention to the low momentum region of validity of the EFT, 
these terms in Eq. (J4.3)) should be negligible compared to those involving fewer derivatives 
of the metric. Their only role is to provide the covariant UV counterterms necessary to 
remove the subleading logarithmic divergences in (Tab). Conversely, if when multiplied by 
finite renormalized parameters a R and j3 R of order unity, they are not negligibly small, then 
the applicability of the EFT framework for low energy gravity is in question. In cosmology 
certainly a necessary condition for this framework to be applicable is that the Riemann 
curvature tensor components and their contractions are negligibly small in Planck length 
units, i.e., 

4, \R ab cd \ = hG N \R ab cd \ « 1 , (4.7) 

and we restrict ourselves to this regime. 

The independent (unrenormalized) components of the energy-momentum tensor with 
non-vanishing expectation values in a general RW initial state are the energy density, 



£« = (Ttt)u 

, 3£ 



4tt 2 



[dk] k 2 <7fc 



y 2 + 



+ m' 



\<t>k? 



2tt 2 



and the trace, 



T 



2tt 2 



[dk}k 2 a k 2#Re(0*0 fc ) + 



[dA;] k 2 <jfc 



(4.8a) 



+ m 2 + £R 



m 



2tt 2 



[dA;] k 2 a k \ 



(4.8b) 



The other non- vanishing components of (T a j,) in a general RW state are the diagonal spatial 
components, (Ty) = pgtj. The isotropic pressure p may be obtained from the energy density 
e and trace T, by p = (e+T)/3. The conservation equation (|4.4jl in the case of RW symmetry 
has only a time component which is non-trivial, namely, 



e + 3H (e + p) = e + H {As + T) = . 



(4.9) 



The unrenormalized expressions (J4.8j) satisfy this relation by use of the equation of mo- 
tion (|2.1ip . provided that a comoving momentum cutoff &m, introduced to render the in- 
tegrals finite, is itself independent of time. An important criterion for any renormalization 
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procedure is that it preserve this property so that ()4.9|) remains valid for the fully renormal- 
ized quantities as well. Notice that a fixed cutoff in the physical momentum pu = kM/a will 
not preserve the form of the covariant conservation Eq. (|4.9|) . because of the non- vanishing 
time derivative operating on the upper limit of the integrals at ku = Pm®, if Pm — M is 
assumed to be independent of time. 

In the case of spatially homogeneous and isotropic RW spacetimes the adiabatic method 
has been shown to be equivalent to a fully covariant treatment of the divergences of the 
energy-momentum tensor which preserves its conservation j^. The starting point of this 
method is the WKB-like form of the exact mode functions, 

M*) = \J^kr k ex p / dt W)) , (4.io) 

which when substituted into 1)2.11)1 . yields the second order equation for Q k , 

^H«-s>-K* + THt-Jt t" 1 ' 

From this expression a systematic asymptotic expansion of the frequency Qk in time deriva- 
tives of the metric scale factor a(t) can be developed. At leading order, neglecting all time 
derivatives, Q k ~ u k . Substituting this into the right hand side of (j4.11)) . one finds to second 
order, 

^ uj k + R-^{H + 3H 2 ) + § ^H 2 + . . . , (4.12) 

where the ellipsis consists of terms third and higher order in derivatives of the metric. It is 
clear that this asymptotic expansion is valid at large k, i.e., at distance scales much shorter 
than the characteristic scale of the variation of the geometry H~ l . Hence requiring the 
exact solutions of the mode equation (J2.ll)) to match this asymptotic expansion to some 
order implies that the quantum state density matrix of the scalar field f)3. lUj) should match 
that of the local vacuum to that order. It is a statement of the Equivalence Principle 
in the low energy EFT that the local, short distance properties of the quantum vacuum 
at a point x should approximate that of the nearly flat space vacuum constructed in a 
local neighborhood of x. Hence the wave functional (JH.22)) must have large k components 
characterized by {( k , ^k} which are universal, corresponding to local geometric invariants at 
x in the effective action, and the same for all physically realizable states, independent of the 
geometry of the spacetime at larger scales. 
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When ()4.12|) is substituted into ([4.10)1 . and the resulting mode function is substituted 
into ()4.8|) with set equal to one, one obtains integrands which match the quartic and 
quadratic divergent behavior of the unrenormalized stress tensor components 38]. Up to 
adiabatic order two these are explicitly given by 

h 



,(2) 



T (2) 



with 



-(2) 



,(2) 



U k + 



m 



H 2 



47r 2 a 3 

h 

47r 2 a 3 



m - k 



[dk] k e 



2„(2) 
k ' 



[dk] k 2 T^ 



(4.13a) 
(4.13b) 
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C H 2 m ff 2 
H — — H 



(4.14a) 



5 m 6 



H 2 



Uk 



H + H 2 + 



m 
2^2 



2H + 3H 2 + 



3 m 

^4 



H 2 



(4.14b) 



Notice that these expressions are state independent and universal, depending only upon the 
RW geometry and the parameters m, £ of the matter Lagrangian. Although not manifestly 



covariant in form, Refs. |47[ show that subtracting these second order asymptotic terms 
from the unrenormalized energy density (|4.8aJ) and trace ()4.8b|) corresponds to adjustment 
of the generally covariant counterterms up to two derivatives in the low energy effective 
action (|4.1|) . Consistent with this, it may be checked that the second order energy density 
of Eq. ()4.14aj) and the second order pressure p^ 2 \ satisfy the covariant conservation 
equation (|4.9|) . provided any cutoff of the k integrals is again independent of time. Hence 
the (partially) renormalized energy density e u — e^ 2 \ and trace T u — T^ 2 \ which are free of 
quartic and quadratic divergences, also obey the conservation equation ([4.9)1 . 

In order to remove the remaining logarithmic divergences in the energy density and trace 
in a general RW spacetime, the terms containing up to four derivatives of the metric must be 
subtracted as well in four spacetime dimensions. The expressions for the adiabatic order four 



terms in the mode expansion QUID , or and can be found in j2y, |4l|. We shall not 
need their explicit form here, and simply assume that one can identify a particular solution 
Vk to (|2.11|) . whose frequency function Qk possesses an asymptotic expansion for large k 
which agrees with ([4.12)1 . up to fourth adiabatic order, and stress tensor components ([4.8)1 . 
which agree with and up to fourth adiabatic order. 
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In the general case, this is the necessary and sufficient condition for the renormalized 

stress tensor to be finite and conserved in the RW state corresponding to this particular 

mode function t> fc . For example, in de Sitter spacetime, these modes, Vk, could be taken to 

n 

be the Bunch-Davies (BD) modes |2j, since these are adiabatic order four modes and the 
BD state is a candidate vacuum state. A general set of modes <pk can be written then as a 
Bogoliubov transformation (|2.19|) of these vacuum modes. The difference of the renormalized 
stress tensor in this general state with that given by the particular choice of <pk = Vk and 
(Tfc = 1 then define the finite state dependent terms in the stress tensor in the general RW 
initial state. In order for the initial state defined by this general set of modes to remain a 
UV allowed RW initial state, the state dependent terms in the renormalized stress tensor 
should not spoil the fourth order approach to the local vacuum which we required of the 
vacuum modes Vk- Hence we impose the condition that the integrals with state dependent 
integrands must be convergent as well. Pure or mixed states satisfying this condition will 
be called UV allowed RW states. 

These UV allowed states are described by mode functions, (pk and corresponding density 
matrix parameters, 717.; <Jk) for which the integrands in the stress tensor components (|4.8jl 
agree with the fourth order adiabatic integrands ejj, and at large k. We may choose 
any particular fourth order adiabatic Vk with respect to which to define the renormalized 
vacuum energy-momentum tensor components, 



-(4) 



<Pfc=ffc 

<Tfc=l 



(4.15a) 



T = T 

- 1 v — u 



- TW. (4.15b) 

<Pk=v k 

<7fc=l 



The definition of the class of UV allowed states then guarantees that the difference of stress 
tensors for any UV allowed RW state with respect to this choice of vacuum are well-defined 
and finite. To identify these terms we have only to introduce the form of the Bogoliubov 
transformation (|2.19J) for the general mode function <pk into 1)4.8)1 . and using (J2.20)) . separate 
off the vacuum terms evaluated at Af. = 1, Bk = 0, and o~k = 1, which are renormalized 
by 1)4.15)1 . The remaining terms are the finite terms for arbitrary UV allowed RW states 
with respect to the given vacuum choice. Collecting these remaining state dependent terms 
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gives the fully renormalized result, 



e = (T tt ) 



2tt 2 



1 



R 



+ 




(4.16a) 



and 



T = (T) 



R 




) 



(4.16b) 



where Nf~ is defined by Eq. ()2.28|) . The vacuum terms denoted by the subscript v defined by 
Eqs. ()4.15|) and the additional state dependent terms in Eqs. f)4.16)) are separately conserved. 
Because the state is assumed to be UV allowed, Nk also approaches zero faster than k~ A as 
k — > oo, and all terms in Eq. (J4.16)) are finite, i.e., there is no cutoff dependence and the 
integrals may be extended to infinity. Note also that a pure state with n& = 0, = 1 remains 
a pure (coherent) state under the Bogoliubov transformation 1)2.19)1 . notwithstanding the 
non-zero value of Nk = \Bk\ 2 for this state in the Vk basis. The quantum coherence effects 
of the Bogoliubov transformation appear also in the rapidly oscillating interference terms 
involving A^B* k in Eq. ()4.16)) . which must be retained in the general UV allowed RW coherent 
state in order to retain the strict time reversibility of the evolution, as we shall see in Sec.lVll 



The development of the previous sections applies to general RW initial states of the scalar 
field of any mass and £ in an arbitrary RW spacetime. In this section we apply this general 
framework to the special case relevant for slow roll inflationary models, namely de Sitter 
space with a massless minimally coupled inflaton field. If spatially flat sections are used, 
then the scale factor for de Sitter space takes the form, 



V. SHORT DISTANCE EFFECTS IN INFLATION 



QdS — 




Hr] 



1 



— OO < 7] < , 



(5.1) 
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with H a spacetime constant related to the scalar curvature by R = 12H 2 . The entire 
de Sitter manifold may be represented as a hyperboloid of revolution embedded in a five 
dimensional flat Minkowski spacetime [32j. The hyperboloid possesses an 0(4, 1) invariance 
group of isometries, which can be made manifest if spatially closed coordinates (e = 1) are 
used. The flat coordinates (e = 0) with the scale factor given by (|5.1|) cover only one half 
of the full de Sitter hyperboloid. None of the results presented in this section will depend 
on the choice of flat, open, or geodesically complete closed spatial sections, so we treat only 
the flat sections (e = 0) in detail. 

In the flat sections under the transformation of variables y = —kr] = kexp(-Ht), the 
mode equation (|2.12|) becomes Bessel's equation with index, 



" 2 = J-^-12«. 



The Bunch-Davies (BD) state 



48, 



49j is the unique RW allowed state which is completely 
invariant under the full 0(4, 1) isometry group of de Sitter space. In the coordinates where 
the scale factor is given by (j5.1J) the BD state is specified by the particular solution of (j2.11j) 
given by 



U D = . / nh c inu/2 c in/ I 



L ~ w \U,r 
1 fh 



H$-\y) (5.3) 



e ik \ as y -> oo . 

a V 2k 

For v 2 < 0, v is pure imaginary and Eq. (|5.4|) is independent of the choice of sign of Im(z/). 
Note that the asymptotic form for y = \krj\ —* oo, holds independently of the value of v. 
From the subleading terms in this asymptotic expansion of the Hankel function for 
large values of its argument, it is straightforward to show that the BD state (with = 0) 
is an adiabatic order four UV allowed RW state for any v. Hence taking Vk = </>f D is 
allowed and the adiabatic order four subtractions of Eqs. (j4.15|) yield a UV finite vacuum 
energy-momentum tensor expectation value in the BD state, which satisfies T v = —4e v or 
p v = —e v , as a consequence of the de Sitter invariance of this state. The calculation of the 
renormalized e v as a function of the parameters m, if, £ is given in Refs. QHQ- 

The massless minimally coupled field is of particular interest both because slow roll 
inflationary models rely on such a field, and because it obeys the same mode equation in a 
RW spacetime as gravitons in a certain gauge (^J. For this field m = £ = 0, v — 3/2 and 
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Eq. ()5.4|) becomes 



Although this state is perfectly UV finite, an infrared divergence occurs in the two-poin ; 
function, Eq. f)2.29j) . The BD state must therefore be modified at very small values of k 51 1 
which means that, strictly speaking, it is not possible to take Vk = 0f D for all k. An IR 
finite vacuum state is the Allen- Folacci (AF) state [28], which is actually a family of IR 
finite states. Because these states are not de Sitter invariant, the energy- momentum tensor 
for the massless minimally coupled scalar field in any of these states is also not de Sitter 



invariant 



28, 



521 ] . Nevertheless it has been proven 27[ that the energy-momentum tensor 



of the m = 0, £ = scalar field for any of the AF states and indeed for any UV finite, 
homogeneous, isotropic state, asymptotically app roaches the de Sitter invariant energy- 
momentum tensor found by Allen and Folacci [2£j, namely, p v = —e v = 119 HH A /9607T 2 . 

Since an AF state is just the BD state modified at very small values of k, the short 
distance or UV properties of the AF and BD states are identical. In this paper we are 
only concerned with these short distance effects, so we take as our preferred vacuum state 
Vk = even in the m = £ = case, ignoring the infrared divergences in the two point 
function which this generates. This is not a problem for the power spectrum provided that 
the k « modes are not the ones that dominate today. There are no infrared divergences 
in the energy- momentum tensor of the BD state and the finite difference in (T a fc)# between 
the BD state and any realistic AF state are of order if 4 and small if H <C Mp\. Hence 
this distinction will play no role in our analysis of short distance modifications of the initial 
state. 

With v k = 0f D the power spectrum for a general UV allowed pure state is given by ()2.32j) 



becomes 



where 



P,(k; t) = + £ (\B k \ 2 \<P B k D \ 2 + M^kBl (0f D ) 2 ]) (5.5) 

"fc=0 



*r=£i*PT=»(i)V*Y). ^ 

is the spectrum of the Bunch-Davies state for the massless, minimally coupled field. In the 
late time limit, rj — > 0~, 
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the BD power spectrum becomes completely independent of k, i.e., scale invariant. If one 
evaluates the power spectrum at the time of horizon crossing instead, k\r]\ = k/aH ~ 1, one 
also obtains a scale invariant spectrum with a normalization differing slightly from (J5.7j) by 



a constant factor of order unity 



20 



The terms in (15. 5j) dependent on \Bk\ 2 and on A^B* k are the contributions to the value 
of the power spectrum for states different from the BD vacuum state. The first important 
point to notice is that if the state is UV allowed then \Bk\ must approach zero and P^ must 
approach P? D at large k. From this fact we can draw an immediate conclusion, namely, 
if Pcf, is evaluated at horizon crossing, k = Ha = e Ht , then the power present in any UV 
allowed initial state always reverts to its scale invariant BD value for fluctuations with large 
enough k which cross the horizon at sufficiently late times. To make this statement more 
quantitative suppose that the initial state is non-adiabatic up to some physical scale M at 
the initial time to with a(to) = do, while above that scale is the same as the BD state. The 
comoving wave number corresponding to this scale is ku = Mao. The horizon crossing time 
for a mode with this wave number is 

t M = H- 1 In k M = t + H- 1 In J . (5.8) 



Fluctuations which leave the horizon at times t > tM will have k > ku and the standard 
BD power spectrum, P^ B . Thus, if inflation goes on for longer than ln(M/H) e-foldings, 
the initial state effects at the physical scale M inflate to scales far outside the horizon. If 
the scales we observe in the CMB now correspond to k > kM, i-£-, to modes which left the 
horizon of the de Sitter epoch at times t > tu-, then there will be no imprint of the short 
distance initial state effects at scale M in the present day CMB observations. Conversely, 
if k < k^f for the currently observable modes then the initial state modifications of the 
spectrum at scale M may be observable. Taking the present horizon crossing scale to be of 
the order of the present Hubble parameter H now , this would imply that the condition 

Ma /a now ~ H now (5.9) 

is satisfied. This condition on Mao/a now is a general constraint on the present observability of 
any initial state effects in inflation in the low energy EFT framework, regardless of their short 
distance origin. Additional constraints and additional parameters may arise in any given 
inflationary model. For example in the slow roll scenario the measured CMB power spectrum 
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depends on the slow roll parameter e in Eq. (|2.34|) . so that observational constraints on the 
CMB power spectrum generally depend on more parameters than simply those of the initial 
state of the scalar field. 

A constraint which does not depend on other parameters of the inflationary model is 
that arising from the energy-momentum tensor of initial states different from the BD state. 
If these contributions to the stress tensor are too large the model will deviate significantly 
from de Sitter space and may not inflate at all. Specializing our general results (j4.16|) to 
the case of de Sitter space with flat spatial sections and a scalar field that is massless and 
minimally coupled, the relevant energy-momentum tensor components are: 

£ = £bd + h + h, (5.10a) 
T = T B d + 2/ 1 -2/ 2 , (5.10b) 

£ + T T /2 

P = = Pbb + h ~ — , (5.10c) 

where the finite state dependent integrals I\ and I2 for the case nu = 0, = 1 are 

h = ^^d^ 2 (|5 fc | 2 |0ri 2 + Re[A fc S*(0f D ) 2 ]) , (5-lla) 
h=^~ 2 /°°dA;A; 4 (|S fc | 2 |0^| 2 + Re[^(0f D ) 2 ]). (5.11b) 

In the following subsections we consider various examples of modifications of the initial 
data for inflation and make use of these general expressions to compute the modified power 
spectrum and energy-momentum tensor components they generate. 



A. de Sitter Invariant a States 



The BD state for the massive scalar field described by the mode functions (|5.4jl is a special 
RW allowed state, since it is invariant not only under spatial translations and rotations, but 
also under the full 0(4, 1) isometry group of globally extended de Sitter spacetime fsj. For 
that reason it has seemed the most natural analog of the Poincare invariant vacuum state of 
QFT in flat Minkowski spacetime, and is usually assumed, explicitly or implicitly, to be the 
relevant "vacuum" state of the scalar field in inflationary models. However, as maximally 
extended de Sitter spacetime is very different from flat spacetime globally, global 0(4, 1) 
invariance is a much stronger condition than local flat space behavior. 
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Because it is a RW allowed state the BD state indeed has a two-point correlation function 
($(a;)$(s')) with short distance properties as x — > x' that depend only on the local geometry 
at x, and more specifically are of the Hadamard form [3^. Since is of mass 

dimension two, this is equivalent to the statement that the BD mode functions (|5.4j) possess 
an asymptotic expansion for large k which agrees with Eq. (|4.1U|) and (|4.12j) up to second 
adiabatic order. However, any UV allowed state satisfies this property. What is special 
about the BD state is that its asymptotic expansion for large k agrees with the adiabatic 
expansion (j4.12J) to all orders. This is a much stronger statement than that it goes over 
to the local Poincare invariant vacuum state in the flat space limit, since an infinite order 
adiabatic state carries information about the geometry of the background spacetime at all 
scales, including correlations on causally disconnected scales much larger than that of the 
horizon if -1 , a situation which has no analog in flat space. 

The fourth order adiabatic condition on the state guarantees that the stress tensor in that 
state possesses no new divergences, and can be renormalized accordingly by the standard 
local counterterms of the low energy EFT of gravity. None of these RW allowed states are 
de Sitter invariant except the BD state. Nevertheless, we showed in Ref. j^f that as a 
consequence of the redshifting of short distance modes to large distances in de Sitter space, 
all RW allowed initial states for Re v < 3/2 have a renormalized (T a b) which approaches 
the BD value at late times. All such fourth order RW states are equivalent locally and are 
a priori equally possible initial states for an inflationary model. 

The only possible way to generalize the BD state while maintaining de Sitter invariance, 



for Re v < 3/2, wou 
and satisfy ipTHty [IE, 



be to require \Ak\ = \A\ and \Bk\ = \B\ be independent of k 
. Because of the unmeasurable overall phase of the Bogoliubov 
coefficients, we can choose A to be purely real and parameterize these squeezed states by a 
single complex number, z = re td as 

1 1 

A k = coshr = — , = = , (5.12a) 

|A| 2 VI - e-+°* V ; 

A e a 

B k = e ld sinh r = — = = , . (5.12b) 

y/l - |A| 2 VI - e a+a * 

The two alternate parameterizations shown in terms of A and e a are sometimes employed 



13.48 



54|. 



The Wightman function (|2.29J) for this general one complex parameter family of squeezed 
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states is de Sitter invariant, in the sense that it is a function only of the 0(4, 1) de Sitter 



invariant distance between the points. However, as pointed out in Ref. |55J] none of these 
squeezed a states are truly invariant under the de Sitter isometry group except for the BD 
state, r = 0, since the r ^ states transform under 0(4, 1) transformations by a non-zero 
(in fact, infinite) phase. Since the Bogoliubov coefficient in (|5.12j) does not approach 
zero at large k, this class of r or a states are not UV allowed states unless r = identically. 
If Eqs. (15.12)1 were taken literally for all k, the integrals ()5.11j) and ()5.10|) would diverge 
quartically. Although these states were used in various contexts, such as studying the sign 
of backreaction effects of partic e creation in de Sitter space Wl, and have been reconsidered 

Finn 

lately by several authors [8J, H3J, |54| , this severe UV divergence is unacceptable for a physical 
initial state within the low energy effective theory of gravity described by Eq. (|4.1|) . This is 
clear even at the level of field theory with no self-interactions, provided only it is covariantly 
coupled to gravity, since the stress tensor in the general a state has divergences which depend 
on a, and thus requires state dependent counterterms for its renormalization 33] . When 
self- interactions are considered still other unphysical features become manifest [54]. 

Computing the power spectrum for a general (r, 9) state by using Eq. (|2.31j) with (|5.12j) . 
we obtain 

P#;r,0)=P* D {l+2smhV-^^ . (5.13a) 

At late times, rj —>■ 0~ (with k fixed) 

hH 2 

PJk; r, 9) -> — - (l + 2 sinh 2 r - sinh 2r cos 9) . (5. 13b) 

An 2 

Because of the non-adiabatic UV modification of the BD state by r at arbitrarily large k, the 
effects of this modification do not redshift away, and essentially the same result is obtained 
if Eq. ()5.13)1 is evaluated at horizon crossing time i] = — 1/k — with krj = — 1 fixed. The 
scale invariant modification (|5.13|) is equivalent to that found in Ref. |8( with the choice, 



sinhr = — — , (5.14a) 



which gives 



2M 
2M 

2 If '~ \2M 
H f2M\ H 2 „ / M\ 1 m ( H (2M 



p * = p r \ 1 - m sin \-h) + U) r p * BD i 1 - F sin [it ) >• (515) 
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where the last approximate equality holds if H <C M, with M the physical UV scale of 
new physics, denoted by A in [$] . In the case of exact de Sitter invariance the general de 
Sitter invariant squeezed state gives a simple multiplicative correction to P^ D sinusoidally 
varying with H/M, which is itself unobservable in the CMB, since it has no k dependence. 
Hence it could be interpreted as a redefinition of the inflation scale H, with no observable 
consequences. 

More importantly, the k independence of the strictly de Sitter invariant r and 9 state 
is untenable in the EFT framework, since it leads to state dependent divergences in the 
energy-momentum tensor. These are avoided if and only if the Bogoliubov coefficient Bf~ 
approaches zero fast enough at large k, for the state to be a UV allowed state. Thus one 
could consider a cutoff version of the (r, 9) states in which r = above some large but finite 
comoving cutoff, 

k M = Ma , (5.16) 

with M the physical cutoff (in units of inverse length) at some arbitrary initial time to. One 
can then assume that modes with k > ku are in the adiabatic BD state while modes with 
smaller values of k are in an (r, 9) state. It is clear that such a state is no longer de Sitter 
invariant and has a power spectrum identical to (|5.13|) or ()5.15|) for k < ku, but reverting 
back to its BD value for k > ku- Thus in such a state, 

= pBD 1 1 + _ ^ ^ ginh 2 r _ ginh 2r cog flj | ^ ( 5A7 j 

instead of ([5. 13)1 . There is now a sharp break in the power spectrum, which could be 
observable in principle, if we are fortunate enough to have access to the right values of k ~ &m 
in the present CMB. If we assume that this condition is satisfied by the wavenumber of the 
present CMB, then observations would put a constraint on the magnitude of the deviations 
from scale invariance of the spectrum of the form (j5.17j) . However, since this constraint is 
model dependent in any given inflationary model, we do not consider it further, and turn 
instead to the constraints arising from the contributions of such a cutoff r state to the 
energy-momentum tensor during the de Sitter phase. 

It is clear from the divergence of the energy- momentum tensor at infinite ku that the 
dominant contribution to the integrals in (|5.11j) comes from those modes close to the UV 
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cutoff. Substituting ()5.12|) into ()5.11j) and cutting the integrals off at ku 
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(5.18b) 



(5.19) 



The properties of the functions Fq (x) are discussed in Appendix El All contributions to 
Eqs. (15.18(1 are finite (for finite M) at all times, including the initial time. All terms redshift 
with the expansion at least as rapidly as a -2 , in accordance with our general theorem in 
Ref. [23| . The terms involving Fq are rapidly oscillating for early times, Mao ^ Ha, but 
redshift to zero as fast or faster than the non-oscillatory terms for late times, Mao Ha. 
The transition from the oscillatory to damping behavior occurs at a time when Mao ~ Ha 
which is of the same order parametrically in M/H as £m, defined in Eq. (|5.8|) . By that time 
all the oscillatory terms give contributions to I\ and I2 which are already of order if 4 and 
negligible. 



Since the maximum of Fg P \x) is of order unity, at an x of order unity, while |_Fg Pj (x)| is 
bounded by 1/x as x — > 00, the oscillating terms are never larger parametrically than HM 3 , 
while the non-oscillating terms make a maximum contribution to the energy density or the 
pressure of order, 

^-sinh 2 r (5.20) 
I67H 

at the initial time, a = ao- Comparing this with the energy density of the inflaton field at 
the onset of inflation, 3/j-if 2 Mp[/87r and requiring that the backreaction from the additional 
terms (|5.18j) be smaller gives the bound, 

H M Pl 



(p)l 



sinhr < v67r 



(5.21) 



M M 

It is possible for the right side of this inequality to be larger than unity, even for H <C M. 
Hence sinh r could be quite large and the break in the power spectrum ()5.17(1 large enough 
to be observable, without creating too large a backreaction. This kind of a non-adiabatic 
initial state modification of the BD state at short distances produces the largest effects in 
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the power spectrum ()5.17|) of the CMB without giving an unacceptably large backreaction 
during inflation. 

If on the other hand r is assumed small, as for example in ([5.14)1 . then no term in ([5.18)1 
is larger in magnitude than 

Ifo^' (5 ' 22) 

and we obtain only the weaker condition, 

M < v^Mp! . (5.23) 

In this case the effects on the power spectrum would be small, though observable in principle, 
and the physical momentum scale of the cutoff ku ja§ = M at the onset of inflation need only 
be somewhat smaller than the Planck scale, beyond which there would be no justification 
for using the low energy effective action for gravity ()4.1[1 in any case. 

Summarizing, the r states do not match the adiabatic expansion of the mode functions 
or the energy-momentum tensor at any k for which r ^ 0. They are therefore completely 
non-adiabatic states. For that very reason there is no bound on the size of their effects on the 
CMB power spectrum for k below the cutoff scale ku- Observations of the CMB may provide 
the strongest constraints on this kind of initial state modification, but the quantitative bound 
depends on the inflationary model. The only model independent constraint for these non- 
adiabatic modifications of the initial state comes from the magnitude of the backreaction 
produced, which is a relatively weak constraint, giving Eq. ()5.21jl or ()5.23jl . The largest power 
of the physical cutoff Mallowed by dimensional analysis appears in the stress tensor, i.e., 
M 4 in ()5.20j) for these states. To illustrate how these results change if adiabatic conditions 
are imposed on the initial state, we consider next zeroth order adiabatic states. 



B. Adiabatic Order Zero States 

A state of given adiabatic order can be obtained by first substituting the expansion ()4.12j) 
into ()4.1()jl and expanding to that adiabatic order. The result, evaluated at some arbitrary 
time to serves as the initial condition for the exact modes (ftk- These modes will remain 
adiabatic to this order for all time |32j ]. 

A zeroth order adiabatic vacuum state for the massless minimally coupled scalar field 
can be obtained by setting Q k = u k in Eq. (J4.1U)) and omitting terms proportional to a in 
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the resulting expression for <fik since they are of first adiabatic order. At the time to one has 
then 

fe (O) = -HrioJ^j-e-**, (5.24a) 

fWk 

Mo) = v k (o) = J y e ~ lkm ■ (5 ' 24b) 

Substituting into ()2.22|) with v k = 0f D gives: 



A k= + (5.25a) 
B t = ^e-^. (5.25b) 
If the relations (|5.25|) are substituted into Eq. (|5.5|) . one finds the late time power spectrum, 



P^-O") =P ( 



BD 



1 _ sin(2/c?7o) sin 2 (fcr? ) 



(5.26) 



In this case the initial state effects produce sinusoidal modulations of the power spectrum 
with wave number which vanish as k — > oo. Also note that if the adiabatic order zero initial 
condition is taken in the infinite past, r] — ► — oo, the modifications vanish as well. This 
is because in that limit the adiabatic order zero initial state becomes the BD state with 
Ak — 1, Bf. — 0, and the spectrum reverts to the standard BD value. 

Danielsson 7] has considered initial data which are of the same form as Eqs. ()5.24j) 
and (|5.25p . but despite this apparently adiabatic construction, rather than viewing r] as 
a fixed time Cauchy surface, where initial conditions are imposed on the state for all k, 
he takes rjo to depend on k in such a way that krjo = —M/H, with M a fixed physical 
scale. The motivation seems to have been to avoid making any statement about modes 
whose physical wavelength is shorter than the cutoff M _1 , and indeed rjo(k) = —M/(Hk) is 
the conformal time at which the mode with comoving wave number k first falls below the 
physical cutoff M. However, inspection of ()5.25|) with this substitution, shows that \B k \ now 
behaves as a constant, H/2M, as k — > oo. Hence this prescription yields a state which is not 
adiabatic at all, but amounts to populating the highest frequency modes considered with 
a constant particle occupation number, and choosing a cutoff (r, 9) state with parameters 
given by (|5.14j) . Thus the results of the previous subsection apply. These initial conditions 
taken literally for all k lead to an energy-momentum tensor which is quartically dependent 
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on the cutoff, just as in the previous subsection. As discussed there, this is not a physically 
allowed UV state if extended to arbitrarily large k, i.e., arbitrarily late times rjo(k) — ► 0~. 

This example illustrates the shortcomings of considering modifications of the initial state 
of inflation and their effects on the CMB power spectrum alone, without also considering the 
associated effects on the energy-momentum tensor and backreaction. When one considers 
only the power spectrum for some finite range of k, it may seem perfectly reasonable to 
restrict attention to only those modes with a physical wavelength larger than the short dis- 
tance cutoff scale M _1 , since no sum or integral over k is required for the power spectrum. 
However, the stress tensor does require such a sum over all k, and some prescription for the 
ultra short distance modes has to be given as these modes will redshift to wavelengths larger 
than the short distance cutoff at later times. The essential question is not at what time 
7/0 these UV modes have physical wavelengths larger than M _1 , but rather what contribu- 
tion do these short distance modes make to the energy-momentum tensor, which involves a 
sum/integration over all k, at any time. This question requires that a choice be made about 
whether the stress tensor is to be a consistent source for Einstein's equations and only if 
it is, can the magnitude of the backreaction be reliably estimated. General covariance of 
semi-classical gravity requires the state to be adiabatic at the very highest trans-Planckian 
energies, and this adiabaticity condition in turn constrains the possible effects of short dis- 
tance initial state modifications on the power spectrum at late times, which might otherwise 
be overlooked. 

Instead of taking t]q to be a function of k let us assume that all the modes are determined 
at the same arbitrary but fixed time rjo, independent of k, by Eq. (|5.24jl . This defines 
a true adiabatic order zero state. Although B k given by ()5.25|) now does decrease with 
increasing k, its magnitude still does not fall off fast enough to make the state fourth order 
adiabatic and UV allowed. Since k 4 \Bk\ ~ k 3 as k — > oo, the energy- momentum tensor 
can depend as much as cubically on the comoving momentum cutoff of the mode sum. The 
cubic divergence in the state dependent mode sum means that there is no local (state- 
independent) counterterm available to absorb this divergence. The necessity of imposing a 
physical cutoff on the behavior of ()5.25|) implies that the power spectrum ()5.26|) cannot be 
valid for arbitrarily large k either, but instead must approach the BD spectrum more rapidly 
than ()5.26|) as k —>■ oo. If we insert a cutoff as in the previous subsection such that the 
modes are the BD modes for k > knt, then the condition (|5.fJj) is necessary for these initial 
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state modifications to be observable in the CMB today. 

If Eqs. (I5.25J) are substituted into Eqs. (|5.11|) with a cutoff ku one finds that 
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As expected the effects of the state dependent terms redshift away like a 2 and a 4 , and are 
largest at or near the initial time rj = rjo when they are of order, 

hHM 3 



32tt 5 



(5.29) 



Requiring this to be less than the energy density of the inflaton field gives the bound, 

M < (UttHM^ , (5.30) 

for the adiabatic order zero state, in place of (|5.21j) for the non-adiabatic state. 

It is clear that Eq. (|5.29|) . softer by one power of H/M compared to the previous 
case (I5.2(J|) is the result of the fact that the adiabatic order zero state has a Bogoliubov 
coefficient, \Bk\, which approaches zero at large k with one power of 1/k in (|5.25|) . If we 
had chosen a state which matches the adiabatic vacuum mode Vk to first, second, or third 
order, i.e., with Bogoliubov coefficient \B^\ approaching zero at large k like k~ 2 , k~ 3 or k~ 4 , 
respectively, then we should expect to obtain leading contributions to the stress tensor com- 
ponents that behave like H 2 M 2 , H 3 M or H 4 ln(M/H), respectively for large M/H. When 
the state is a UV allowed state, the stress tensor components are independent of the upper 
limit ku = MaQ of the mode integrals for large M, so that the integral may be extended to 
infinity. In that case the stress tensor components are of order if 4 , independent of the cutoff 
M, and negligible compared to the energy density driving the inflation for all H Mpi. 
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C. Boundary action approach 



The authors of Refs. [25j] have discussed setting conditions of the form, 

(d t ^k + ^k)\ t=to = 0, (5.31) 

on the initial state mode functions, motivated by the addition of boundary terms to the low 
energy EFT action functional. Here k is in general a complex function of k. For the BD 

state, 

" ~ TTWo (5 - 32) 

becomes purely imaginary in the limit k\r]o\ — > oo. If we make use of (|2.22j) we find that the 
Bogoliubov coefficients are given by 

Ak ~ 2v /(Im«)(Im«BD) 1^(0)1' { } 

B k = (5.33b) 
2V(Im«) (Im« BD ) 1^(0)1' 1 ^ 

up to an overall phase. 

If attention is restricted to modifications of the BD state corresponding to the lowest 
dimension local operator in the scalar EFT on the initial time surface at t = tn, namely 
/5(Vj$) 2 /M where M is again the physical cutoff scale, then the authors of Refs. [25j] argue 
that this would lead to a modified initial condition of the form (|5.31|) with 

(3k 2 

« = «bd + -ytj ■ (5-34) 
%M 

If the effective action on the boundary is real for real time it would seem that (3 must be 



real. The authors of Refs. 



25j | treat (3 as a real parameter, obtaining corrections to the (real) 
power spectrum which are linear in (3. Treating (3 as an arbitrary complex parameter, we 
obtain the Bogoliubov coefficients for the case k\rjo\ ^> 1, 

- -^mj- ik lml3 ) " • (5 - 35a) 

I^-S^^ 11 "")" 1 ' (5 - 35b) 



As in the a state case discussed in Sec. IV Al these Bogoliubov coefficients are non-adiabatic 
and would lead to a divergent stress tensor if continued to arbitrarily large k. Substituting 
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the cutoff ku = Ma® into the power spectrum (|5.5|) as before, we obtain at late times, 

P,(k; „ - 0") = i? D {l + l^^lj (W/Te 2 ^) + \(3\ 2 A cos 2 ^)) } . (5.36) 

If \/3\ is of order one then both terms in Eq. 1)5.36)1 are of the same order at k ~ /cm, and the 
(large) deviations from a scale invariant spectrum may be observable in present CMB data. 
The same remarks about fine tuning to k ~ ku and dependence on the specific features of 
the inflationary model apply to this initial state modifications as to the cutoff r states of 
Sec. IVAl 

If |/3 1 is assumed to be much less than unity, we can write j3 = \/3\e n and obtain the 
modified power spectrum to linear order in \/3\, 
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The finite state dependent contributions to the stress tensor given by the integrals ()5.11)) 
are easily written down for the case of general complex j3, but because of the denominators 
in ()5.35)) they are rather complicated. In the case that (3 is purely real the integrals simplify 
and may be evaluated in terms of the functions = Fq P ^ with the latter defined in 
Eq. (JOl. The result is 
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All terms are again finite at all times and redshift at late times at least as rapidly as a~ 2 , 
in accordance with our general theorem in Ref. 2?j]. The oscillatory integrals F^(x) are 
discussed in Appendix B, and the illustrative particular case of F^ 4 \x) is plotted in Fig. ^ 
Since the maximum of F^ p \x) is of order 3 to 4 at x = kuij] — Vo) ~ 1 (for p even) or 
x = (for p odd), the maximum value of either the terms linear or quadratic in (3 is of order, 
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FIG. 1: The oscillatory function for p = 4 defined by Eq. (|B5|) . as a function of x = fejv# (r/— 

in contrast to (|5.18jl . where the oscillatory terms were smaller than the non-oscillatory ones. 
If the maximum of the contributions (|5.<ffl|) are required to be smaller than the energy density 
driving the expansion, and 8 <C 1 then the terms linear in f3 give the largest contribution 
and the strongest bound, viz., 



If 8 < 1 this bound on the term linear in 8 is of the same order of magnitude as that 
obtained in Re, Q, but d 18 a gr ee 8 with the bound in Re, Q, whose authors Mg ue that 
the term linear in /3 gives no bound on the bulk stress tensor away from the initial boundary 
surface at t — t . 

The autho, of Re, Q equate inters 8 „ch a 8 ^ with a Gaussian cutoff iu t, tathe, 
than with a hard cutoff at k = used here. However this by itself does not account for the 
disagreement. A Gaussian cutoff amounts to a slightly different choice of initial state which 
is also perfectly UV allowed, since it approaches the BD state faster than any power of k. 
This different state gives a different finite contribution to the renormalized stress tensor, 
which is of the same order in M as the contribution (|5.39jl . differing only in its numerical 




(5.40) 
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coefficient at the initial time, and with a smoother, less strongly oscillatory behavior in time 
than that shown in Fig. ^for the hard cutoff hyi- in neither case does the contribution to 
the energy-momentum tensor fall off exponentially in the conformal time difference t] — r] 
as claimed in Ref. [nj]. At late times rj — ► 0~ either the Gaussian or hard cutoff of the 
initial state momentum integral yields a finite stress tensor with components that fall off as 
a -4 , a -3 and a~ 2 , as expected by the redshift of the RW expansion and in accordance with 
Ref. {27! . Hence the finite terms in the renormalized stress tensor which are first order in 8 
are no more subject to renormalization ambiguities or localized on the boundary than those 
which are second order in 8. Which terms give the stronger bound on 8 depends entirely 
on the values assumed for the parameters kfttWo, H and M. 

For general values of 8, if the ratios H/M and M/Mp\ are the same order of magnitude, 
then requiring that the maximum contributions (|5.39jl be smaller than the energy density 
driving the expansion yields the bound, 

max(/3,/3 2 ) < 10 «-> 100 . (5.41) 

If (3 < 1 this bound is easily satisfied, while if 8 > 1 the maximum value comes from the 



terms quadratic in 8, and the results of Ref. 17j are recovered. We conclude that although 



the precise evaluation of the 
estimates given in either or 



backreaction effects of initial state modifications differs from the 
the qualitative final conclusion that the backreaction 
constraints are not a very severe restriction on the parameter (s) of the boundary value action 
is similar to the conclusions reached by these authors. 

VI. ADIABATIC PARTICLE CREATION AND DEPHASING 

One of the principal physical effects that can be described by the semi-classical Gaussian 
density matrix is particle creation by a time varying RW scale factor (5^ . Although the def- 
inition of a "particle" is intrinsically non-unique in a time varying background, it is possible 
to use the adiabatic nature of the UV allowed RW states to constrain this non-uniqueness 
considerably. The adiabatic particle concept was studied in some detail in Ref. j^J, where a 
proposal was made for the adiabatic basis. The deficiency with that earlier proposal is that 
the two parts of the stress-energy tensor corresponding to vacuum and particle contributions 
are not separately conserved. Here we remedy that defect and in the process remove almost 
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all of the non-uniqueness in the definition of adiabatic particle number. 

Let us remark first that it is always possible to express the density matrix p in the time 
independent number basis in which the number operator a^a^ of Eq. (|2.23|) is diagonal. The 
transformation to this particle number basis may be derived by the methods of Refs. 0, 57 1 
with the result, 



where n labels the set of integers {n^}, one for each distinct k. In this occupation num- 
ber representation, the density matrix is time independent and diagonal, since the o~\~ are 
constants of motion. The positive diagonal matrix elements of p in this discrete number 
representation may be viewed as the probabilities of finding exactly particles in the mode 
labelled by wave number k, with the particle number basis defined by the time independent 
operator, = a^a^, i.e., 



which is equivalent to Eq. ()2.23|) . together with ()3.6a|) . 

Since the unitary transformation to the basis exactly undoes the action of the time 
evolution operator, the preceding definition of particle number is always time independent, 
no matter how rapidly the geometry changes with time. Hence it carries no information 
about particle creation in the time evolving RW geometry, or indeed about any features 
of the time evolution of the system whatsoever. Moreover, if one makes a Bogoliubov 
transformation from one exact set of eigenfunctions <pk{t) to another exact set Vk{t), then 
the particle number with respect to the new basis remains exactly time independent, 
cf. , (J2.28|) . Useful as these are in the description of the density matrix and energy-momentum 
tensor of the field, any such time independent parameterization of particle number is not 
the appropriate one to describe particle creation in time varying RW cosmologies. For that 
purpose one needs to define an appropriate time dependent number operator with respect 
to an approximate adiabatic vacuum basis, which agrees with the exact number basis only 
at very large k. 

The definition of the approximate adiabatic basis is constrained by several physical and 
practical requirements. First the basis should be defined by a set of WKB-like approximate 
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mode functions, 

k = a-*f k = ^^4^ ex P (-* / dt'W k {t')-iS k (t)\ , (6.3) 

similar in form to Eq. ([4.10)1 . for some adiabatic frequency function W k and additional time 
dependent phase 6 k (t), to be specified. If W k = fife, with fl k a solution to Eq. (|4.11j) . and 6k 
is constant in time, then f k would be an exact solution to the mode equation ([2.12)1 . and we 
would obtain the time independent number basis (jbMj) . However, there is no constraint on 
the (ftk to satisfy the equation of motion (j2.11|) . except asymptotically for large k, where they 
should approach the exact mode functions. Hence we have considerable freedom to choose 
the two functions, W k (£) and 6k{t) which define the adiabatic basis. The most physically 
meaningful choice of basis is that which requires the vacuum energy-momentum tensor 
defined by <p k to agree with the adiabatic expansion of the vacuum terms in the covariant 
energy-momentum tensor to a fixed adiabatic order. It is this tensor which couples to 
gravity and is covariantly conserved (as distinguished from the canonical Hamiltonian), so 
that matching the adiabatic basis (|6.3|) to this tensor will guarantee separate conservation 
of the vacuum polarization and particle contributions to the stress-energy. 

It is clear that the leading (adiabatic order zero) terms in Eqs. (|4.14j) are vacuum po- 
larization terms, and not particle contributions, since they appear even in static or flat 
spacetimes, where there can be no particle creation whatsoever. Thus, the frequency W k 
of the adiabatic modes ()6.3)) should certainly match the zeroth order term u k in order for 
the energy-momentum tensor of these modes to match the adiabatic order zero terms in 
Eq. (|4.14j) . However, adiabatic order zero matching is not sufficient. If we do not require 
matching <p k to at least second order in the adiabatic expansion of the vacuum v k , then we 
will find particle contributions to the energy-momentum tensor which diverge quadratically 
in the UV cutoff kw These divergent adiabatic order two terms are clearly part of the 
state independent vacuum polarization contribution to (Tab), an d not the particle content 
of the state. Choosing the adiabatic frequency and phase in ()6.3|) to match the adiabatic 
order two expansion terms ()4.14j) exactly guarantees that they will be completely removed 
by the subtractions in (|4.15j) . It is this physical requirement that the power law divergences 
in the conserved stress tensor components should be associated with the state independent 
geometrical contributions and not the state dependent particle content that determines the 
adiabatic basis functions <fr k and renders the definition of adiabatic particle number (almost) 
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unique. The only remaining non-uniqueness of the particle definition arises from the possi- 
bility of matching the vacuum contributions to the stress tensor to higher than second order 
in the adiabatic expansion. 

In order to define the adiabatic basis precisely let the exact mode functions be expressed 
in terms of the adiabatic modes (|6.3j) as 

fk = a k (t)f k + (3 k (t)f*, (6.4a) 

in terms of time varying Bogoliubov coefficients, a k (t) and f3 k (t). The time dependent 
Bogoliubov coefficients are required to satisfy 

K(t)| 2 -|/3 fc (t)| 2 = l, (6.4b) 

for each k, in order to guarantee that the transformation of bases is a canonical one. As 
in the case of the time independent change of bases discussed in Section HU this allows a 
two real parameter freedom in the choice of the f k for each k (up to an overall irrelevant 
phase). Rather than using the Bogoliubov coefficients or W k and 5 k as the two independent 
parameters, it is more convenient to use W k and a different independent function of time 
V k (t) , defined by the relation on the first time derivative of Eq. ()6.4a|) , namely, 

f k = (~iW k + a k (t)f k + (iW k + yj P k (t)ft . (6.4c) 

The canonical transformation from the exact mode functions f k to the approximate adiabatic 
functions f k is now completely specified by W k and V k , while 5 k of Eq. ()6.3|) is fixed (implic- 
itly) in terms of these two. Before determining W k and V k explicitly from the stress tensor 
components we first compute the particle number Af k and density matrix in the general 
adiabatic basis determined by Eqs. ()6.4j) . 

If the original field operator $(t,x) is expanded in terms of the approximate adiabatic 
mode functions, f k and f£, the corresponding annihilation and creation operators, 

a k (t) = a k a k (t) + a\[3* k (t) , (6.5a) 
al(t) = 4<(t) + a_Mt) , (6.5b) 

are generally time dependent. Because of Eq. (J6.4b|) and the freedom to choose the phase 
of a k we may set 

ot k {t) = coshr fc (t), (6.6a) 
(3 k (t) = e^W sinhr fc (t), (6.6b) 
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in analogy with Eq. ([2.260 . With respect to this time dependent adiabatic basis, the number 
density of particles in mode k is 

M k {t) = (a[a k ) = n k + a k \f3 k (t)\ 2 = sinh 2 r fc (t) + n fe cosh2r fc (t) , (6.7) 

which also depends on time in general. However, if W k and V k are chosen properly at large k, 
A4 will be an adiabatic invariant with respect to the effective classical Hamiltonian ([A10J1 . 
and therefore slowly varying in time for a slowly varying RW scale factor. The number 
density of particles spontaneously created from the vacuum with n k = is \(3 k (t)\ 2 . An 
expression for this quantity in terms of W k and V k is easily obtained from Eq. ([6.4(1 in the 
same way as ([2.22(1 . with the result, 



\h(t)\ S 



sinh 2 r k {t) 
1 



iw k n k 



2hW k 



fk + iW k 



2 , 1 / T/ , ^fc 




(6.8) 



where the last expression follows from inserting the WKB-form ([4.10(1 for the exact mode 
function f k in terms of the exact frequency Q k . From this we observe that if (and only if) 
the adiabatic frequency matches the exact frequency, W k = Q k and V k = —Cl k /Q k , then the 
Bogoliubov coefficient f3 k vanishes identically. In that case there is no particle creation and 
A4 = n k is strictly a constant of the motion. 

Using the value of \f3 k \ 2 the density matrix may be expressed in the adiabatic particle 



basis, by the methods of Refs. 



44. 



571 ] . Unlike Eq. ([6.1(1 the off-diagonal matrix e 



44, 



ements of p 



57J. Although 



do not vanish in this basis, and are quite complicated in the general case 
the diagonal elements are time dependent, they depend on time only through the function 
r k (t). This means that they are relatively much more slowly varying than the corresponding 
phase variables 6 k (t), upon which the off-diagonal elements of p also depend. Thus, in this 
adiabatic number basis it becomes possible to argue that particle creation is related to phase 
decoherence or dephasing of the state: since macroscopic observables are generally relatively 
insensitive to the process of averaging over the rapidly varying phases, one can replace the 

|. In the pure 



exact p in this basis by its more slowly varying diagonal elements only 
state case, with adiabatic vacuum conditions in the infinite past, these diagonal elements 
simplify and are given explicitly by 



Pn h =2l k {k]t) 



(2l k 



1)!! 



sech r k tanh 2Zfc r k . 



2'* y 



(6.9) 
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The diagonal matrix elements are non-vanishing only for even integers, corresponding to 
the fact that particles are created from the vacuum in pairs. The positive numbers ()6.9j) 
have the interpretation of the probabilities of finding exactly l k pairs of adiabatic particles 
at time t in the mode labelled by k = |k|, with vacuum initial conditions as t —>■ — oo. If this 
replacement of the exact pure state density matrix (j3.1U|) by its phase averaged diagonal 
elements is justified, then the von Neumann entropy suffers the replacement, 

5 = -Tr(plnp)^-— / dA; k 2 ^ p 2lk (k; t) \np 2lk (k; t) , (6.10) 

71 i k =o 

which becomes time dependent. Although, in general the effective von Neumann entropy 
does not grow strictly monotonically in time, starting in an initial pure state with all of 
the rfc = leads to a larger effective entropy at late times when some of the modes have 
rk ^ 0,0,0]. This shows that particle creation is directly related to increased squeezing 
of the initial state, and the growth of entropy this entails corresponds to the effective loss of 
information resulting from averaging over the rapidly varying phases e ±t6k< ^ in macroscopic 
physical observables. 

The validity of a truncation of the density matrix to its diagonal terms only in the 
adiabatic number basis and the associated loss of phase information will depend on the 
initial state and the details of the evolution [58]. When particle creation takes place from 
initial vacuum-like states this would seem to be quite a good approximation in the cases 



that it has been tested quantitatively 
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44j ]. Conversely, if one starts from a different 



type of initial state the squeezing parameters need not increase monotonically with time, 
and phase averaging is not justified. If it should happen in some special case(s) that the 
squeezing coefficients are constant in time for all k, so that the adiabatic particle number 
basis becomes an exact vacuum basis, then by making the appropriate time independent 
Bogoliubov transformation (j2.19|) to that basis one can set all the = 0. Then it is clear 
that the 9 k become undefined and no phase decoherence of the initial state can occur by 
particle creation or dephasing effects in 9k- 

With these general remarks on the adiabatic particle number basis and dephasing let us 



fix the still undetermined functions Wk and Vk- Following Ref. |26J] let us replace the exact 
mode functions by (pk, and their time derivatives (pk by d(pk/dt = (—iWk + Vk/2 — 3H/2)(f)k 
with a k = 1, p k = 0, and n k = in Eqs. ()4.8|) for the stress tensor components. Using ()6.4|) 
we obtain the (cutoff dependent) adiabatic vacuum contributions in this basis to the energy 
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density and trace, namely, 



h 



T 



47r 2 a 3 
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with e k and T k the following functions of W k and V k 



[dk] k e k , 
[dk] k 2 f k 



(6.11a) 
(6.11b) 
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(6.12a) 
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The adiabatic vacuum basis and particle number may be fixed now by setting these two 
expressions in terms of W k and V k precisely equal to the corresponding state-independent 
vacuum contributions to energy-momentum tensor in a general RW spacetime to second 
adiabatic order, given previously by Eq. (|4.14|) . That is, we require 



fk 



.(2) 



T 



(2) 



(6.13a) 
(6.13b) 



thus giving two relations for the two functions, W k and V k which determine the adiabatic 
particle basis. Solving for W k we find 



m 2 + (6^ - 1) [2uo 2 + HV k + H)+(6£- l) 2 [HV k + H+% 



2(6^-l)4 2) +T^ 

If this is substituted back into either of Eqs. (|6.13j) . we obtain a quadratic equation for V k , 
so that the two functions W k and V k which determine the adiabatic vacuum modes can be 
determined algebraically for any £ and m for the general RW geometry. 

In this way the energy-momentum tensor components may be written in a form analo- 
gous to (|4.16jl . where because of the exact matching ()6.13|) the vacuum terms up to second 
adiabatic order are now identically zero. Therefore the remaining non-vacuum terms de- 
fined with respect to the time dependent adiabatic vacuum necessarily satisfy the covariant 
conservation equation (|4.9|) . The non- vacuum terms in the adiabatic particle basis can be 
obtained from ()4.16|) with the replacements of N k — > J\f k , v k — > <f) k , A k a k and B k — > (3 k . 



,(2) 



(6.14) 
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The non-vacuum energy-momentum tensor components for general Wk and Vk are given 
explicitly in j^, and we do not to repeat them here. 

Since Eq. ()6.13|) matches the vacuum energy-momentum tensor to second adiabatic order, 
the adiabatic frequency Wk obtained from solving Eqs. (J6.13)) agrees with the adiabatic 
expansion (|4.12j) up to and including second adiabatic order, differing from only at 
adiabatic order four. On the other hand since V k 2 and HVk appear in Eqs. ()6.13|h and Vk 
contains terms of odd adiabatic orders, V& need agree with — Cl^/flk — —i^k/^k only up to 
first adiabatic order, i.e., 



general non-vanishing at adiabatic order three, and the lowest order term in the adiabatic 
expansion which appears in the expression for \/3k\ 2 in Eq. ()6.8|) is sixth order, i.e., the adia- 
batic particle number J\f k defined by Eqs. f|6.7|) . ()6.8|) and (j6.13J) is a fourth order adiabatic 
invariant. Time derivatives of A4 are correspondingly highly suppressed, and particle cre- 
ation is small in slowly varying backgrounds, particularly at the highest wave numbers. In 
flat spacetime, H, e and Vk vanish, and (jQ4j) together with (j4.14j) give Wk = ojk which is 
time independent. Hence the adiabatic modes become exact modes, [3k = 0, and no particles 
at all are created. 

In the general case, Eqs. ()6.13)) match the vacuum energy-momentum tensor contributions 
to second adiabatic order and therefore satisfy a weaker condition than the fourth order UV 
allowed state condition. However the second order condition (j6.13|) is sufficient to render 
the total number of adiabatic particles created finite, i.e., 



Note that this would not be the case had we matched the energy-momentum tensor com- 
ponents only to the lowest (zeroth ) or der adiabatic expansion, or followed the Hamiltonian 
diagonalization procedure of Ref. [56]. The requirement of matching the stress tensor con- 
tributions 1)6.12)1 exactly to a fixed adiabatic order as in ()6.13)1 (rather than re-expanding 
the algebraic expressions for Wk and Vk) removes essentially all of the ambiguity in the 
definition of the particle number, and guarantees that the non-vacuum particle contribution 
is separately conserved. 




(6.15) 



where the ellipsis includes terms of third adiabatic order and higher. Hence Vk + flk/^k is in 




(6.16) 
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Because of the weighting of the momentum integrals for the energy-momentum tensor 
components (j4.16|) by an extra power of k with respect to Eq. (j6.16|) . the identification (j6.13|) 
still leaves a logarithmic cutoff dependence in the vacuum energy-momentum tensor com- 
ponents. Unlike the power law divergences which are non-covariant in form, the remaining 
logarithmic cutoff dependence is proportional to the geometric tensor of Eq. (|4.6a|) . 

Hence it can be viewed as a geometric vacuum polarization contribution in the low energy 
EFT equipped with a short distance cutoff. Even if the short distance cutoff is of the order 
of the Planck length, the higher derivatives in mean that this contribution is negligible 

for all geometries varying more slowly than the Planck scale. Of course one is free to define 
the adiabatic particle basis by replacing the second order terms in (|6.13|) with the corre- 
sponding adiabatic order four (or higher) components and which would remove 
also the covariant logarithmic cutoff dependence in (T a b). Matching ()6.3|) to fourth order 
would make the resulting particle number ()6.7|) an adiabatic invariant of even higher (in fact, 
eighth) order. Matching to higher orders requires accurate knowledge of higher and higher 
order time derivatives of the metric scale factor, and hence of the entire spacetime evolution. 
Because only the power law divergences are non-covariant in form in the adiabatic proce- 
dure, matching to higher than second order is both unnecessary from the point of view of the 
covariant UV divergence structure of the stress tensor, and contrary to an EFT approach in 
terms of the minimal number of spacetime derivatives whenever EFT and the limits ()4.7|) 
apply. The tensor ^H a t actually vanishes in the important special case of de Sitter space, 
which renders the distinction between matching to second or fourth adiabatic order to re- 
move the remaining logarithmic divergences in the stress tensor superfluous in this case. 
Finally, the second order definition of particle number through Eqs. ()6.7|) . ()6.8|) and (j6.13J) is 
also the minimal one necessary to yield a finite total number of created particles f)6.16j) in a 
general RW spacetime. We emphasize that such a definition of particle number is intrinsic, 
based on the physical requirement of identifying and removing the vacuum contributions to 
the conserved stress tensor, and does not depend critically on the existence of flat in and 
out regions of the spacetime, or any extraneous notion of particle detectors Q]. 

We now apply this general second order definition of adiabatic particle number to two 
important special cases with zero mass. When m = the relations (|6.13|) simplify consid- 
erably. For the massless conformally coupled field (m = 0, £ = 1/6), the condition on the 
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trace (|6.13b|) becomes empty but the first condition (|6.13ajl gives directly, 

m = 0, (W k -u k ) 2 + {Vk ~ H)2 = 0, (6.17) 

6 4 

which is solved uniquely for real W k and V k by 

W fc = w fc =-, (6.18a) 
a 

Because the fourth order adiabatic terms in and T^ 4 ^ vanish for the massless conformally 
coupled field in an arbitrary RW spacetime, this result for W k and V k remains unchanged 
if the fourth order adiabatic basis is used. Indeed as a result of the relations (|6.18|) . the 
adiabatic basis functions ()6.3|) become exact solutions of Eq. (j2.11|) for the massless, confor- 
mally coupled field in a general RW spacetime. This means that the Bogoliubov coefficients 
a k and f3 k become time independent. As in Eq. (I2.19J) . M k may be identified with the time 
independent N k and there is no production of massless, conformally coupled scalar particles. 

In a second important special case, the massless minimally coupled field (m = 0,£ = 0), 
V k drops out of Eq. (|6.14|) . and either of Eqs. (|6.13|) gives a quadratic equation for V k , which 
is easily solved. Thus for m = and £ = we obtain 

2(k 2 -e) k 

W k = J- , 6.19a 

2k 2 - e + a?{H + 2H 2 ) a' v ' 



V k = 3H- 



2 ( k% e) ~ 2 \m 2 k 2 -(a 2 H + e)(H + 2H 2 -^ 



(6.19b) 



2k 2 -e + a 2 (H + 2H 2 ) 

in a general RW spacetime. Although this definition of the adiabatic basis was determined 
by matching only to second adiabatic order, by Eqs. ()6.13j) . in the special case of de Sitter 
space it becomes exact. Explicitly in flat spatial sections, e = 0, for which H = 0, we have 

k 2 k 

m = 0,£ = 0,a = a dS : W k = , (6.20a) 

K + ti a dg ads 

_ k 2 + 3H 2 a 2 s _ W k 
Vk ~ H k 2 + H 2 a 2 s -~W k > (6 - 2 ° b) 

with a^s given by Eq. (|5.1|) . As a result of these relations, the adiabatic mode function (p k 
is an exact solution of the mode equation ()2.11j) in de Sitter space. Indeed from Eqs. (|5.1|) 
and fEU} , 

j\t'W k {t') = k 3 J" ^£ a =k V - t*a-\kn) ■ (6.21) 

57 



Hence, 

' W 



2a 3 W h 



exp (-i J &t'W k (t')^j 



TT h(l + k 2 rf) ( ., . ,„ . 
V 2k? eX ^ v % ~^ % n 



m=0 
«=0 



(6.22) 



by Eq. (|5.4|l . The fact that the second order adiabatic mode functions <pk of Eq. ()6.3j) coincide 
with the exact BD mode functions up to a constant phase for the massless, minimally 
coupled field in de Sitter space implies that the Bogoliubov coefficients a k and /3& are time 
independent. This may be verified explicitly by computing Eq. (J6.8|) with W k and given 
by Eqs. ()6.2()j) . Thus A4 may be identified with N k , and there is no production of massless, 
minimally coupled scalar particles in the special case of exact de Sitter spacetime. This 
implies that there is also no phase decoherence of the free massless inflaton field in the de 
Sitter epoch, with respect to the adiabatic particle basis defined by Eqs. (I6.13J) . In the 
next section we will corroborate the absence of decoherence for the massless inflaton by 
computing the decoherence functional directly. A comparison of this result with the results 
of earlier work, such as that of Ref. is given in Appendix C. 

To conclude this section we remark that the adiabatic particle number basis should 
provide an efficient approximation of the low energy semi-classical limit of the energy- 
momentum tensor. If we neglect the phase correlated bilinears, 

(a k a k ) = <TfcSinhr fc cosh r k e~ tdk , (6.23) 

which should make a relatively small contribution to (T a &) in the mode sum over k, compared 
to the terms involving Af k , then the energy density becomes simply 



J [dk]ee ( k 2) Af k . (6.24) 



e ~ e = e v + 

Z7T 

This should provide a useful analytic approximation to the energy density in a general UV 
allowed RW state whenever the RW scale factor varies slowly enough for the higher order 
adiabatic corrections to and A4, to be negligible. 

(2,') 

The quantity e ( k ' in Eq. JESP, 

given by (|4.14aj) has the interpretation of the single particle 
energy in the time varying RW background. Because this second order single particle energy 
satisfies 

4 2) = -#4 2) -#7f , (6.25) 
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where T k is given by Eq. ([4. 14b)) . and the finite vacuum contributions of Eqs. (|4.15jl are 
separately conserved, it follows that the conservation Eq. (|4.9|) is exactly satisfied, provided 
that the trace in the quasi-particle approximation corresponding to Eq. ()6.24|) is 

T^f = T v + 1 ±- 3 J [6k] k 2 if Af k - J [dk ] P 4 2 ) M k . (6.26) 

Defining the ideal fluid pressure in this approximation to be that in the absence of any 
particle creation, i.e., 

P - ^ + ^3 / ^ k2 (4 2) + Ti 2) )N h , (6.27) 
the conservation Eq. (|4.9|) becomes then 



i + 3H (e + p) = 




(6.28) 



after transposing the last term of Eq. ()6.26|1 to the right hand side of Eq. (16.28)) . Thus in the 
quasi-particle limit where it is valid to make the replacements (16.24)) and (|6.26j) . the term 
involving Mk carries the interpretation of the rate of heat dissipation per unit volume due to 
the non- conservation of adiabatic particle number Mk- If the particles are in quasi-stationary 
local thermodynamic equilibrium at the slowly varying effective temperature T e $(t), then 
this rate of heat dissipation may be equated to T e g times the rate of effective entropy density 
generation s e g, 

r -TT = aP? /l d *] fc M 2 '^. (6.29) 
by the first law of thermodynamics. The effective entropy generation gives rise to an effective 
bulk viscosity in the energy-momentum tensor due to particle creation, even in the absence 
of self-interactions of the quasi-particles. 

Let us emphasize that the entropy generation and bulk viscosity are only effective, i.e., 
an approximation valid only to the extent that the phase information contained in the off- 
diagonal elements of the exact density matrix (|3.10|) in the adiabatic particle basis cannot be 
recovered. Likewise the correlations ()6.23j) . which also depend on the rapidly varying phases 
Ok conjugate to Mk should make a negligible contribution to macroscopic physical quantities. 
If the exact phase information is retained, then the evolution remains unitary, as required by 
the equivalence of the evolution to that of the effective classical Hamiltonian (jA10)) . However, 
the extension of the usual description of the cosmological fluid by non-ideal terms is suggested 
by the adiabatic particle creation rate and phase averaging in the low energy EFT. This may 
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provide a useful phenomenological description in some circumstances, and also shows the 
approximations which are necessary in principle to pass from the fully reversible field theory 
description of matter in cosmological spacetimes to an effective, irreversible kinetic theory 
with a definite arrow of time. 



VII. DECOHERENCE 

The density matrix description of the evolution of arbitrary RW states in Sec. IIHI al- 
lows us to describe the quantum to classical transition, i.e., decoherence, in a cosmological 
setting. The fundamental quantity of interest is the decoherence functional between two 

- n 

different histories, Ti 2 or Ti 2 |35(, given by Eqs. (|3.26|) or ()3.27|) . of the pure or mixed state 
cases, respectively. Evaluating the Gaussian integrals needed to compute Tr(pip2) with the 
measure (|3.18|) gives 

1 f un ,2, J 1 ( (2 k , Clk\ 2 . (Clfc^fc — (2k^lk) 2 



12 



^ (<7lfc — (T2k) ■ {7A) 



!r7 2 ./ 1 J llV C,/.. G.kJ !> 2 

4 

where {Cifc, ^ik] &ik} and {(2k, ^2k\ o"2fc} are the Gaussian state parameters of the two dif- 
ferent histories. This simplifies somewhat in the case of two different pure state histories, 



Imri2 = r i2 



l =cr fe2 = 1 



' f [dk] k 2 In { -(^- + Y + (Clfc7F2fc 7 ^*^*)' I . (7.2) 
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In this pure state case the decoherence functional may be expressed in terms of one complex 
frequency function, 

T k = A " 7^ ( 7 - 3a ) 

= -ia 3 % - 6i£Ha 3 , (7.3b) 
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where the last relation is derived in Appendix A. It is straightforward to rewrite Eq. ()7.2|) 
in the two alternate forms, 
1 



ImT 
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[dfc] r In 



IT I T* |2 
|Jj_fc_+ j_2fcl 

4(ReT lfc ) (ReT 2/ 



(7.4a) 
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[dife] A; 2 ln{ \a\ 



J 2k 



(7.4b) 



Because of the infinite product of integrations in the functional measure (|3.18|) . it is clear 
that some condition(s) will have to be imposed on the two states or density matrices in these 



expressions, in order to insure a convergent result for large k. 



34, 



he otherwise ill-defined di- 
The divergences in 



vergent nature of ()3.26|) has been noted by several authors 
the decoherence functional are similar to those encountered in the unrenormalized expres- 
sions for the energy- momentum tensor components, Eqs. (|4.8|) . In that case, the superficial 
degree of divergence was quartic, whereas that of Eq. (|7.1|) or (|7.2|) is reduced by one power 
of k, and can be no more than cubically divergent at large k. A method to handle this 
cubic and subleading linear cutoff dependence of the decoherence functional 1)7.1)) or ()7.4)) is 
needed before meaningful results in the low energy EFT description can be obtained. 

In order to study the cutoff dependence of the inner product and decoherence functional, 
let us write 



Tifc — Tfc + 5Tk , 

2k = i-k — I k , 

and expand the logarithm in Eq. ()7.4a)) to second order in STk- We find that 
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(7.6) 



The leading behavior of Re at large k may be read from the first order equation satisfied 
by this function, viz., 

e 



u 2 k + (6t-l)(6ZH A + 



m£Ha 6 T k . 



(7.7) 



which is easily derived from the definition (J7.3a)) and the relations (j3.30cj) and (j3.30d|) . From 
Eq. (|7.7J) we see that 
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Since 5T k is the same order as generically at large k, the decoherence functional ()7.6|) 
will generally diverge as the UV cutoff is removed. Indeed the asymptotic behavior given 
by Eq. (|7.8|) implies at large k that 

5T k 2ka (5a) + ... (7.9) 

so that ()7.6|) will diverge cubically in the cutoff kM, unless 5a = identically. Hence no 
meaningful comparison between two different RW geometries in the low energy EFT can be 
made through the decoherence functional. 

This divergent short distance behavior of the decoherence functional could have been 
anticipated from the relationship between T 12 and the Closed Time Path (CTP) 37J action 
functional. Variation of the CTP action functional with respect to the metric, g a b, pro- 
duces connected correlation functions of the energy-momentum tensor with particular time 
orderings of their arguments. The first variation is the same as that of S eS [g] in Eq. ifO]) . 
and produces the stress tensor expectation value (|4.H|) . The second order variation (|7.fi|) is 
formally proportional to 

Im f d 4 x f clV ^5g ab (x) Tl abcd (x, x') 5g cd (x') ■ (7-10) 

Jo Jo 

We shall be particularly interested in the specific homogeneous variation, 5g a b = 2(5a/a)g ab 
which preserves the RW symmetries. The symmetrized expectation value of stress tensors, 

lmU abcd (x.x') = - (T ab (x)T cd (x') + T cd (x')T ab (x)) (7.11) 

I — I 

is proportional to the noise kernel of fluctuations around from the mean (T ab ) [60]. In flat 
spacetime, where one may transform conveniently to the momentum representation (|7.10j) 
is proportional to the cut in the one-loop polarization diagram of Fig. corresponding to 
the squared matrix element for the creation of particle pairs by the perturbation 5g ab . 

In coordinate space Eq. ()7.11|) is a singular distribution at coincident points x = x', 
involving in general up to four derivatives of 5^'(x,x'). Thus, the second variation ()7.10|) 
exists only for metric variations, 5g ab , which fall off rapidly enough in both space and time 
to permit integration by parts of the derivatives of 5^(x,x'). Since the time interval [0, t) 
is finite in Eq. f)7.1()j) . the surface terms generated by these integrations by parts do not 
vanish at the endpoints in general, and can generate up to two derivatives of delta functions 
at equal spacetime arguments, i.e., formal cubic (and subleading linear) divergences in the 
decoherence functional, which is just what is obtained in from Eq. (|7.6|h (j7.8J) . and ()7.9|) . 

62 




FIG. 2: The imaginary part of the one-loop vacuum polarization given by Eqs. (|7.10|) and (|7.11|) . 
which enters the decoherence functional at second order in the metric variation, 5g a b (represented 
by the wavy lines). The shaded part of the diagram represents the complex conjugation of the 
unshaded part, and the resulting squared amplitude is proportional to the probability for the 
creation of a scalar particle/anti-particle pair, represented by the diagonal cut through the diagram. 

The unrenormalized decoherence functional is not only cutoff dependent in general but 
also ambiguous in that it depends upon the precise definition of the local measure 
which does not affect the physical content of the evolution described by the density ma- 
trix In fact, had we used the = aq^ field to parameterize the density matrix rather 
than gk, then in order to normalize the state properly the measure would have to be replaced 
by a product of the dq^. This differs from Eq. f!3.18|) by an infinite number of local factors of 
a{rf). Hence, we should expect the inner product and decoherence functional denned by this 
measure on the field configuration space to differ from the previous one by cutoff dependent 
contact terms. Indeed, the frequency function would be modified to 

T* = ^|, (7.12) 

cr 

for which the leading term in the large k limit (|7.9)1 cancels. The explicit form of 5Tk is 
given below by Eq. (|7.22|) . with the result that the decoherence functional ()7.4aj) generally 
diverges cubically as the comoving momentum cutoff ku —>■ oo. Because of the different large 
k behavior of if the same steps leading to Eq. (J7.22j) are carried out for instead, the 
resulting expression lacks the last term in Eq. ([7.22)1 . anc 



for that reason its contribution to 
the decoherence functional is only linearly divergent [3^, 3(|. This shows that the degree of 
divergence is dependent on the field parameterization and the definition of the inner product 
(13.26)1 on the field configuration space, which should not have any physical consequences for 
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decoherence due to a slowly evolving geometry in the low energy EFT. 

Because of the appearance of divergences of odd powers, associated with the boundaries 
of the region of integration, the divergences in the decoherence functional cannot be removed 
by renormalization of the bulk terms in the low energy effective action (J4.1)) . Instead their 
renormalization requires introducing counterterms in the effective action, of dimension one 
and three, which are strictly localized on the boundaries. Adiabatic regularization may be 
used to define the necessary subtractions of the decoherence functional, corresponding to 
renormalization of these boundary terms, in a way quite analogous to the adiabatic sub- 
tractions used to define the renormalized energy-momentum tensor jssj . Since terms in the 
effective action up to dimension three are involved, we define the renormalized decoherence 
functional by subtracting from its unrenormalized value the adiabatic expansion of the func- 
tional up to and including its third adiabatic order asymptotic expansion for large k. Since 
the decoherence functional involves an absolute square, this requires subtracting only up to 
the first adiabatic order in the expression inside the absolute value signs in Eq. ()7.4|) . That 
is, we define the renormalized decoherence functional by 
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(7.13) 



where the subscript 1 denotes the expansion of the quantity in parentheses up to and in- 
cluding first order in its adiabatic expansion. This leaves the leading unsubtracted behavior 
at large k to be second adiabatic order under the absolute value signs, and fourth adiabatic 
order in its square in T 12 . This corresponds to subtracting all surface divergences up to and 
including third adiabatic order in the CTP action functional, which is what is required. 

To carry this out explicitly we begin by rewriting Eq. (j7.3b|) in conformal time in the 
form, 
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and computing its first variation, 
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The variations of Xk and its derivative are computed by varying Eq. ()2.14|) . to obtain 
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which is solved in terms of the retarded Green's function of the differential operator on the 
left hand side, 
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(7.17) 



in the form, 
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(7.19a) 
(7.19b) 



The lower limits of the integrals in Eq. f)7. 19j) depend on the initial conditions of the wave 
functional and give a time independent phase in the decoherence functional below which 
we shall not need to specify. Substituting Eq. ()7.18|) into ()7.15j) and using the Wronskian 
condition (J2.17)) . we find that the 5a k term cancels and we are left with 
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Using Eq. f)2. 17|) again, we have 
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ReT fc Xk V ^ 
This form is valid for a scalar field of arbitrary mass and curvature coupling in a general 
RW spacetime. In order to renormalize it, we should subtract its asymptotic expansion up 
to adiabatic order one and substitute the square of the subtracted quantity in Eq. (|7.13jl . 
We carry out this subtraction explicitly in two important special cases, namely, when the 
mass m = and the curvature coupling £ takes on either its conformal or minimally coupled 
value, £ = 1/6,0, respectively. 

Using Eq. (|7.22j) together with ()7.19|) and the form of Eq. ()2.14j) in the massless, confor- 
mally coupled case, we obtain 
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which involves no time derivatives and hence is clearly of adiabatic order zero. If we were to 
substitute this directly into Eq. ()7.4a|) . we would obtain a cubically divergent decoherence 
functional. Since this cubic divergence can be removed by simply redefining the canonical 
variables as in Eqs. (jA16|) and (|7.12j) . it is clear that it can have no physical significance. 
However, if we first subtract off the adiabatic order zero part (which in this case is the entire 
expression), then the renormalized decoherence functional (j7.13J) in an arbitrary RW space- 
time is identically zero for the massless, conformally coupled field. This lack of decoherence 
corresponds to the lack of particle creation for this field in any RW space which we have 
found in the previous section. 

In the massless, minimally coupled case the corresponding expression is 

_ f v A , x fa"\ r , Al2 2i x fa'\, |2 , 2i X * k 
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which is generally non-zero. Note that the first term in this last expression depends upon 
the variation of the RW scale factor over the entire evolution from an arbitrary (unspecified) 
initial time at the lower limit of the 7/ integral to the final time 77, while the last two terms 
of Eq. (|7.24j) depend upon the variation of the scale factor only at the final time. It is 
these two latter surface terms that generate divergences in the unsubtracted decoherence 
functional (JZHJ), when integrated over k. Such surface terms arise in the covariant expres- 
sion (|7.1U|) if the conservation of T ab is used to express the tensorial noise correlator (|7.11|) 
in terms of covariant derivatives of scalar quantities, and then an integration by parts is 
performed. 

Specializing to de Sitter spacetime and using the form of the BD mode functions (|5.4|h 
the last term in Eq. ()7.24j) becomes 

2i X k 
ftXk 

Here the first term is of adiabatic order zero as in the previous conformally coupled case, 
while the second term is of adiabatic order one. Hence both terms are fully subtracted in the 
renormalized decoherence functional. The other contact term becomes in de Sitter space, 
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which consists of an adiabatic order one and order three term. Hence subtracting up to 

adiabatic order one removes the first term but leaves the l/k 2 rf term unaffected. Since the 
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first term in Eq. (|7.24|) involves two time derivatives of the scale factor or its variation, it is 
adiabatic order two, and likewise unaffected by the subtraction of up to adiabatic order one 
terms. Hence finally, 
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which is UV finite and non-zero for arbitrary variations of the scale factor. 

If we consider the particular variation of the scale factor in which the de Sitter Hubble 
parameter H is varied in Eq. (|7.27|) . while the conformal time rj is held fixed, then a' /a = 
— 1/f] and a" /a = 2/rf are fixed and Eq. ()7.27j) vanishes. Hence we find that under variations 
of the de Sitter curvature, which compare the wave functionals of the quantum field in 
macroscopically different de Sitter universes but at the same conformal time, the decoherence 
functional for the massless, minimally coupled field vanishes identically. 

The ambiguous contact terms which are field parameterization dependent are removed 
by the adiabatic regularization and renormalization procedure in Eq. ( j7.13jl . as the two 
cases considered explicitly above show. Hence, Eq. (|7.13j) yields both a finite decoherence 
functional free of unphysical dependence on the short distance cutoff, and one that is inde- 
pendent of redefinitions of the scalar field variables and inner product. The renormalized 
decoherence functional proposed here vanishes in the two special massless cases of the con- 
formally coupled field in a general RW background and the minimally coupled field in a de 
Sitter background, the same two cases studied in the previous section where the adiabatic 
particle creation rate Mk = 0. Since the imaginary part of the polarization tensor (j7.11j) 
is just the cut one-loop diagram shown in Fig. which is proportional to the probability 
for the metric fluctuation to create a particle/anti-particle pair from the vacuum, a close 
correspondence between the lack of particle creation and a vanishing decoherence functional 
for homogeneous metric perturbations in the Hubble parameter is not unexpected. The fact 
that the adiabatic subtraction procedure for the decoherence functional proposed here sup- 
ports this correspondence suggests that it is the correct one to define a finite decoherence 
functional for semi-classical cosmology. In order to prove that this is indeed the unique 
procedure for defining a physical decoherence functional in the EFT approach, the adiabatic 
subtractions of first and third order should be related to definite boundary counterterms 
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in the CTP effective action which reside exclusively on the surfaces at the initial and final 
;imes. These boundary terms may be related to those found recently by the authors of Ref. 
I9I ] . We leave the determination of these surface terms for a future investigation. 



VIII. SUMMARY AND CONCLUSIONS 

The principal objective of this paper has been to place semi-classical cosmology within 
a consistent EFT framework, in which possible short distance effects can be parameterized 
by well-defined initial conditions at the onset of inflation. Although the general principles 
and adiabatic methods underlying such an EFT framework have been available for some 
time, we have thought it worthwhile to make these assumptions completely explicit in this 
paper, and demonstrate how they can be applied and extended in a number of different 
ways, which may be useful for future cosmological models and observations. Because of 
the several different applications considered in the paper, we collect here and summarize 
the principal results, together with the relevant equations and sections where each point is 
discussed in detail. 

• The general homogeneous, isotropic RW pure state is defined by field amplitudes (pk 
obeying ()2.11|) . which are linear combinations of vacuum modes ()2.19|) with Bogoliubov 
coefficients satisfying ()2.2()|) . 

• These pure RW states are squeezed vacuum states annihilated by ak in the mode ex- 
pansion ()2.6|) and specified by two real time-independent squeezing parameters (j2.26J) . 
up to an overall irrelevant phase. 

• The wave functionals of these pure RW states are Gaussians ()3.22|) in the Schrodinger 
picture field coordinate basis. 

• The general RW state with a non-zero occupation number ()2.23|) is a mixed state 
described by the mixed state Gaussian density matrix (j3.10J) in the coordinate basis. 

• The general RW mixed state requires three independent functions of k, (Ck,^k',o'k), 
which are related to the mode function 0^ by ()3.6|) or ()3.7|) . and determine the three 
equal time symmetrized correlators of the field by ()3.8|) in the Hamiltonian description. 
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The first two of these functions of k, 71"/-) are time dependent and form a canonically 
conjugate pair for the unitary evolution of the density matrix ([3.29)1 . ()3.30)) . described 
by the effective classical Hamiltonian ([A10J) - ([A12|) . in which ft appears as a parameter. 

The third function of k, at = 2^^ + 1 is strictly a constant of the motion. 

The form of the Hamiltonian of the scalar field evolution in a fixed RW background 
depends on the field parameterization, and in general is not equal to the covariant 
energy density e = T tt . 

The covariant and Hamiltonian descriptions of the evolution are completely equivalent 
nonetheless, and the total Hamiltonian of the combined matter plus geometry system 
(|3.34j) vanishes by time reparameterization invariance for evolutions satisfying the 
classical Friedman equation. 

The power spectrum of scalar field fluctuations in the general homogeneous, isotropic, 
mixed RW state is given by Eqs. ()2.32|) . 

The spectrum of the actual scalar metric fluctuations observed in the CMB are de- 
pendent upon additional parameters which may be different for different inflationary 
models. An example is the dependence on the slow roll parameter e in Eq. (|2.34|1 . 

The energy density and trace of the stress tensor in the general RW state is given by 
Eqs. (14.16)1 . with e v and T v the values in the fiducial vacuum state. 

In order to be a UV allowed RW state with short distance behavior consistent with 
general covariance of the low energy EFT and the Equivalence Principle, the fiducial 
vacuum state and all other physical states must be fourth order adiabatic states. 

Any modification of the fourth order adiabaticity condition at short distances has the 
potential to disturb the conservation of (T a b), and/or violate the Equivalence Principle 
at arbitrarily large distances and late times, which would also violate the decoupling 
hypothesis of low energy EFT. 

The Bunch-Davies (BD) state is a UV allowed fourth order adiabatic state, which is 
also invariant under the full 0(4, 1) isometry group of global de Sitter spacetime. 
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The general one complex parameter squeezed a states of the scalar field in de Sitter 
space are not UV allowed fourth order adiabatic states (even for non-self-interacting 
scalar fields), except for the single value of the parameter corresponding to the BD 
state. 

Because all UV allowed states are fourth order adiabatic, their power spectrum ap- 
proaches that of the BD state for sufficiently large comoving wavenumbers k > ku-, 
and sufficiently late times t > after the onset of inflation ([5.8)1 . when EFT methods 
should apply. 

As a consequence of this kinematic effect of the expansion, any modifications of the 
power spectrum due to initial state effects require a coincidence of fine tuning ([5.9)1 in 
order to be observable in the CMB today. 

Assuming such fine tuning and cutting off the squeezed a state at a finite large comov- 
ing momentum scale ku produces potentially observable scale dependent modifications 
of the CMB power spectrum (|5.17|) . whose magnitude depends in general upon addi- 
tional parameters of the inflationary model. 

Cutoff a states and non-adiabatic states generally produce the largest backreaction 
contributions during the inflationary epoch, given by ([5.18)1 . which are of order ([5.20)1 . 

States which are adiabatic order zero up to the cutoff scale ku produce scale depen- 
dent modifications of the CMB power spectrum ([5.26)1 which may be observable as 
modulations in the CMB power spectrum. 

Such states also produce backreaction effects during inflation which are somewhat 
smaller in amplitude than the cutoff a states, and which can be calculated exactly 
from Eqs. (EHUll . lETTTl) and (l5~2Tll . 

The modifications of the initial state given by the addition of a local higher dimension 
operator with coefficient (3 in the boundary action approach are non-adiabatic and 
yield in general modifications to the CMB power spectrum at linear order, Eqs. ()5.36j) 
and 1ET37I) in (3. 
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The backreaction contributions to the stress tensor during inflation are also linear in 
(3 in general and of order /3M A , which may be significant, depending on the cutoff and 
inflation scales M and H, but do not disturb inflation if (|5.4U|) is satisfied. 

The adiabatic expansion of the stress tensor can be used also to define a time dependent 
particle number basis for particles created by the RW expansion, (|6.7|) . with parameters 
Wk and Vk defined by Eqs. fl4.14j) . (jdl2j) and (jfi.lHj) matched to the stress tensor 
exactly at second adiabatic order. 

The total particle number defined in this way is the minimal one that is finite, (|6.16jh 
and gives separately conserved vacuum and particle contributions to the covariant 
stress tensor. 

In the general massive case the particle number is not conserved but is a sixth order 
adiabatic invariant, implying that the density matrix in the adiabatic particle rep- 
resentation has slowly varying diagonal components but much more rapidly varying 
off-diagonal components. 

Although the exact evolution is unitary and fully reversible, averaging over the rapidly 
varying off-diagonal elements of the density matrix in this basis gives rise to an ap- 
proximation which is effectively dissipative, and in which the effective von Neumann 
entropy of the reduced density matrix ()6.10|) may increase with time. 

Neglect of these same phase correlations in the energy-momentum tensor via the ap- 
proximations ()6.24|) and (|6.27|) gives an effective rate of heat dissipation due to particle 
creation ()6.28|) and ()6.29|) . even in the absence of matter self- interactions. 

Notable exceptions to this dephasing occur in several special massless cases, due to 
the absence of particle creation for a conformally invariant scalar field in any RW 
spacetime, and for a massless, minimally coupled scalar field in de Sitter spacetime. 

The latter result implies that the quantum phase information in the density perturba- 
tions derived in slow roll inflationary models is not washed out by the expansion alone, 
so that the loss of phase decoherence in such models must be due to other effects not 
considered in this paper. 
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• The decoherence functional for arbitrary mixed Gaussian states given by Eq. (j7.1|) is 
related to the noise kernel, or imaginary part of the second variation CTP effective 
action flTTUJ) . 

• The adiabatic method may be used again to define the renormalized decoherence 
functional f)7. 13j) in semi-classical cosmology, which is independent of the short distance 
cutoff and field reparameterizations. 

• This renormalized decoherence functional vanishes in the special cases where there 
is no adiabatic particle creation, corroborating the close connection between particle 
creation, dephasing and decoherence. 

• Comparison of the decoherence functional defined here with a previous result in the 
massless case is given in Appendix C. 

• Verifying this adiabatic subtraction through a covariant subtraction of the surface 
terms in the CTP action functional would allow the study of decoherence effects and 
the quantum to classical transition quantitatively and reliably in general RW cosmolo- 
gies for the first time. 
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APPENDIX A: EQUATIONS FOR (k AND n k AND THEIR EFFECTIVE 

HAMILTONIAN 

In this appendix we compute the three independent and symmetric Gaussian variances 
($ 2 ), ($11^ + n^$) and and derive the equations of motion for the density matrix 

parameters ( k and 7Tfc defined in Sec. IHII 

The square of the defining relation for in Eq. ()3.6|) is 



Hence using Eqs. ()2.9|) and ()2.29|) we find directly for the first variance at coincident space- 
time points, 



which is explicitly real and independent of x. In order to compute the second variance we 
differentiate Eq. (jAljl to obtain 



which is the second of relations (|3.7jl . Hence, using Eq. ()3.3|) the second symmetric variance 
at coincident points is 
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By squaring Eq. ()A3j) and using the Wronskian condition (|2.17jl in 
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we obtain 
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which is the third member of Eq. (|3.7j) . Hence the third Gaussian variance is 



„6 



2tt 2 



dA; A; (jfc 



2tt 2 



dkk' 



b k \ 2 + 12£#Re (0 fc 0£) + 36^10 



2tt 2 



dA; k 2 + 



(A7) 



This establishes Eqs. (|3.8|) of Sec. 11111 

The second order differential equation for ( k may be derived by differentiating ()A3jl . and 
making use of Eqs. (|A6j) and ()2.11|) to obtain 



k 2 - 



m 2 + ^ c 



(A8) 



which is Eq. ()3.32|) of the text. 

The equations for the parameters of the Gaussian density matrix may be compared with 
those arising from the purely classical Hamiltonian ()3.4|) . viz., 



n $ = 6^n $ -a 3 



" A3 + e +m 2 + (6^-l) /f 
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er v cr 

Note in particular that for h ^ the equation of motion for w k (j3.30d|) differs from Eq. (jA9b|) 
of the purely classical evolution by the last centrifugal barrier-like term in Eq. (J3.30dj) 
which is a result of the uncertainty principle being enforced on the initial data through the 
Wronskian condition ()2.17|) . 

The first order evolution equations for the parameters ((f), p; Cfci^fci a k) ma y a l so be re- 
garded as Hamilton's equations for the effective classical Hamiltonian, 
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is the effective Hamiltonian describing the Gaussian fluctuations around the mean field for 
the Fourier mode k. It is straightforward then to verify that Hamilton's equations for this 
effective classical Hamiltonian (in which h appears as a parameter), viz., 



dp 



(A13a) 



* = -lfc> (A13b) 

4 = ^ , (M3c) 

ir k = , (A13d) 

are identical with Eqs. ()3.30|) of the text. Hence ( k and n k are conjugate variables with 
respect to the effective classical Hamiltonian (jAlOjl . 

If we define the complex frequency T k by (|7.3a|) of the text, then by differentiating that 
definition and using Eq. (J3.3U)) we obtain its equation of motion (J7.7|) . On the other hand 
Eq. (|A3J) with a k = 1, together with the Wronskian condition (|2.17j) imply 

h - lie? Cdk = -2m 3 Ml , (A14) 

so that dividing by 2£| = 2<p k (j)* k , and using the definition of 7i k in Eq. (|3.6|) we obtain 

Tk = A _ = _ m 3 k _ Qi^Ha 3 , (A15) 



which establishes Eq. ()7.3b|) of the text. 

Finally we remark that 7iq> depends on the choice of variables used to represent the scalar 
field. Indeed, if we choose the conformal field variable, 

X = a$ , (A16a) 

instead of $ and define the conjugate field momentum, 

Tl x = ^ = x'--X (m = 0,e = 0) (Al6b) 

for the massless, minimally coupled field, the canonical Hamiltonian, 

n x = V + ^U 2 + - x n x = an* + d$n $ , (m = o, e = o) , (Ai7) 

z z a 

differs from 7i$ defined by Eq. ([3.4)1 . and neither is equal to the time component of the 
covariant stress tensor which couples to gravity for general m and £. This is to be expected 
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since the canonical transformation from ($, n$) to (x, n x ) is a time dependent transfor- 
mation, and neither Hamiltonian is a conserved quantity. This has the consequence that 
while every representation is physically equivalent, describing exactly the same physical time 
evolution, there is no spacetime or field coordinate independent meaning to the basis which 
diagonalizes the instantaneous canonical Hamiltonian in a particular set of coordinates [56] , 
and no reason to prefer any such choice over any other as a physical particle basis. 



APPENDIX B: EVALUATION OF INTEGRALS 

In evaluating the integrals I\ and I2 in Eqs. (J5.11|) which contribute to the energy density 
and pressure of the massless, minimally coupled scalar field, we encounter integrals of the 
form, 

/km t,2n+l 
dkk 2n sm(2ku -9) = (-yJ^L-F^ikMu) , (Bla) 

J dk k 2n+1 cos(2ku -9) = (-) n ^f TT F^ n+l \k M u) . (Bib) 

Integrals of this kind are easily evaluated by successive differentiation with respect to u of 
the elementary integral, 

kM • m cos6> - cos(2k M u - 9) , „ „ * 

dfc sin(2fcw - 9) = ^—^ = k M F e (k M u) , (B2) 



where 



Thus, in Eq. (IB1J) we have 
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and, in the particular case = 0, 



We also define 



F^\x)= I dyF(y)= I dy S -^ . (B6) 

Jo Jo y 

For any 9 and p > 0, Fq P \x) are damped oscillatory functions whose maxima occur at x = 
for p odd and on the first oscillation for p even. The values of x and F^ at the maximum are 
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of order unity. The leading behavior as x — > oo is obtained by differentiating the oscillatory 
numerator only, i.e., 

F ° P){X) ^ SinX Sin(x " 6) + ° (i) • (B7) 

Thus the absolute value of Fg P \x) is bounded by l/sc for x ^> 1. Hence the maximum 
value of the integrals (jBl|) are of order k 2 ^ +1 and /c^ +2 respectively for u ~ /c^ 1 , while they 
behave like k]% and respectively, multiplied by a rapidly oscillating function of /cm«, 

as kjyfU — > oo. The form of F^ 4 ^(x) as a function of x in the special case of p = 4 and # = 
is given by Fig. [T] of the text. 



APPENDIX C: COMPARISON OF SQUEEZING IN DIFFERENT BASES 

In Sec. IVII we found that in the second order adiabatic particle basis defined by Eqs. 
(Ifi.lHJ) there is no large squeezing and no particle creation of a massless, minimally coupled 
scalar field in exact de Sitter spacetime. In Sec. I VI II we corroborated the lack of true 
decoherence for this field. In this Appendix we compare this result to earlier work, in 
particular to Ref. 4(|, whose authors use a field amplitude "pointer basis." 

Let us first consider a zeroth order adiabatic basis rather than the second order adiabatic 
basis defined by Eqs. (jfi.lHj) . The zeroth order basis function is obtained by replacing W k of 
Eq. (|fj.2()jl by oj^ = fc/aas, resulting in 

r(o) = / h ( ./"*,., , \ I / h 



exp f-i J dt'uj k (t') \ = i y ^j- exp (-ikrj) (CI) 



instead of Eq. ()6.22|) . When the exact BD mode is expressed as a linear combination of 
these zeroth order adiabatic modes in the form of the Bogoliubov transformation, 

rf^aP^+ff, (C2) 
we now obtain the non-trivial time dependent coefficients, 

"I" = (l - ^) , (C3a) 

^ = —^(-2,^). (C3b) 

Clearly in this adiabatic zero basis sinnr^ = = l/\2kr)\ which does not fall off fast 

enough at large k to be fourth order adiabatic, and which also approaches infinity as r\ — > CT. 
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The problem with the zeroth order adiabatic basis (jCl|) is the "particle" number (j6.16|) . 



2w2jik e\^ = —J ik ^oo (04) 

is divergent at any finite time r\. Corresponding to this linear divergence in the total num- 
ber, the energy-momentum tensor of these "particles" is quadratically divergent at large k. 
Clearly this quadratic divergence is a residual divergence of the vacuum stress tensor and 
has nothing at all to do with physical particles, which are unambiguously well-defined in 
the UV limit k — > oo. In that limit of wavelength much smaller than the horizon scale, the 
region of spacetime the modes sample may be approximated as flat, and the effects of the 
time variation of the geometry should be negligible. The rj~ 2 factor in Eq. (jC4|) shows that 
this mismatch only grows more severe at late times, as the zeroth order basis (|C1|) becomes 
more and more different from the second order basis (jfi.22)) . 

The quadratic vacuum divergences can be subtracted only by matching the energy- 
momentum tensor to second adiabatic order as in Eq. (jfi.lHJl . When that is done the particle 
number is finite and no squeezing, phase decoherence or particle creation effect at all is 
obtained. Moreover, since the second order adiabatic modes ()6.22|) are already exact for 
the massless, minimally coupled scalar field in de Sitter space, any ambiguity in the particle 
concept at wavelengths of the order of the horizon scale is irrelevant here. Going to higher 
orders in the adiabatic expansion will not change the result obtained at second order. 

The authors of Ref. 46] define a field "pointer basis" for the mass ess, minimally coupled 
field in de Sitter space. Comparing Eqs. (15), (19), and (47) of Ref. 4(| with relations (jC3|) 
above, we find that the Bogoliubov coefficients, and (3k of Eq. (15) in Ref. 4(| are precisely 
equal to and respectively, of the zeroth order adiabatic basis given by (jC3j) . As we 
have seen, the squeezing parameter with respect to this basis, sinh rjj. = l/\2kri\ — > oo does 
become very large for superhorizon modes in the late time limit. The authors of Ref. 
argue that the large squeezing in this basis leads to an effective decoherence of modes of the 
scalar field much larger than the de Sitter horizon. 

It is clear that the squeezing is very much dependent on the basis in which it is computed. 
In the second order basis ()6.22jl . determined by the structure of the short distance expan- 
sion of the covariantly conserved stress tensor through Eqs. (|6.13jl . there is no mixing of 
positive and negative frequency modes, and no large squeezing of superhorizon modes. True 
decoherence of these modes should occur through other effects, such as those considered in 
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6l| at the time these modes reenter the horizon. 
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